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Abstract 

We identify and resum corrections associated with the kinematic recoil of the hard 
scattering against soft-gluon emission in single-particle inclusive cross sections. The method 
avoids double counting and conserves the flow of partonic energy. It reproduces threshold 
resummation for high-py single-particle cross sections, when recoil is neglected, and Qt- 
resummation at low Qt, when higher-order threshold logarithms are suppressed. We exhibit 
explicit resummed cross sections, accurate to next-to-leading logarithm, for electroweak 
annihilation and prompt photon inclusive cross sections. 
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1. Introduction 



A large class of hard-scattering cross sections in QCD are factorized into convolutions of parton 
distributions and fragmentation functions with hard-scattering functions [|IJ. Important and 
representative examples are Higgs production and Drell-Yan cross sections, at measured invariant 
mass Q and transverse momentum Qt- We shall refer to these reactions collectively as electroweak 
annihilation. At fixed, large Qt, electroweak annihilation cross sections are written in collinear- 
factorized form as 



in terms of evolved, nonperturbative distributions (densities) (pa/Ai^,^) of parton a in hadron A, 
and hard-scattering functions a a ^y (Q 2 /s, Q^/Q 2 , Q 2 //x 2 , /i, a s (/i)), computed as power series 
in a s . Here s = x a x b S is the partonic invariant mass squared, while fi is a factorization scale, 
which for the time being we equate to the renormalization scale. General single-particle inclusive 
(1PI) cross sections for photons and light hadrons at high pt take a similar form, including a 
fragmentation function. 

Many hard-scattering functions have been computed to next-to- leading order (NLO) in a s (fj,). 
Analytic calculations of still higher-order contributions [|| to a are as yet too complex to carry 
out, except for fully inclusive processes, such as the Drell-Yan production of lepton pairs at 
measured invariant mass ||. Nevertheless, general arguments show that the functions a are 
infrared safe to all orders 

Starting at NLO, the computation of a involves cancellations between soft gluon emission and 
virtual corrections. These cancellations produce plus distributions and delta functions, which 
require integration against smooth functions, such as parton densities. The finite integrals, in 
turn, are potential sources of numerically large corrections at each order in perturbation theory. 
Because of their connection to soft-gluon emission, however, such corrections can sometimes be 
resummed to all orders in perturbation theory. 

For example, in Eq. (ffl), a ab ^v includes distributions that are singular at partonic threshold, 
s = Q 2 , where partons a and b have just enough invariant mass to produce the observed final state. 
Defining z = Q 2 /■§, we find at nth. order singularities as strong as a" [(1 — z)^ 1 ln 2n_1 (l — z)] + . 
Threshold resummations, which organize these distributions, have been developed for a large 
class of cross sections [§, ||, |7j. Although these singularities are manifest in cr ab ^ v , they do not 
generally result in large logarithms in the physical cross section, because they are smoothed by 
the integrals over x a and x b in Eq. ([J). Thus, threshold resummation is not a summation of 
kinematic logarithms in the physical cross section. It is rather an attempt to quantify the effect 
on the physical cross section of a well-defined set of corrections in a to all orders. 

Threshold singularities are not the only singular distributions encountered in the computation 
of a. In addition, the perturbative cross section is singular up to a™[(l/Q T ) ln 2n_1 (Q T / Q 2 )]+ in 
a, Eq. (JJ), when the transverse momentum, Qt, of the electroweak boson is small compared to 
its mass, Q || [|. At each order, Qt is balanced by soft gluons, and singularities in the differ- 
ential cross section at Qt = reflect collinear divergences in a not eliminated by factorization. 
These divergences, resummed or not, cancel in the QT-fntegrated cross section, even before the 
integrals over the partonic fractions x a and x b , although the remainder is still singular at partonic 
threshold. 

Much of the recent interest in soft-gluon recoil effects has centered on the normalization and 



dQ 2 dQ 2 T 



a = / dXa <t>a/A(Xa, ft) / dx b (fib/si^b, /■*) CTab^V 
T ab J 



(Q 2 /s,Q 2 T /Q 2 ,Q 2 /fi 2 ,^a s (fi)) , (1) 
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Pr-dependence of single-particle inclusive cross sections |L0| |, particularly direct-photon produc- 
tion at fixed-target energies [fTl] , |12| [13] , 14|, |15 , 16|, |i~7H . The formalism of Qt resummation for Eq. 
([!]) is not immediately applicable to inclusive high-p^ cross sections, because in this case most of 
the transverse momentum of the observed particle is recoil against other high-p^ particles, while 
only a small portion is from soft radiation. A rough-and-ready approach to soft-gluon radiation 
is to introduce intrinsic transverse momentum for the partons in factorized expressions like Eq. 
([!]), typically in the form of an energy-dependent Gaussian smearing of standard parton densities, 
which enhances the cross section. This method, however, certainly involves double counting, and 
does not respect the conservation of partonic energy. Some time ago, Li and Lai explored the 
possibility that nonperturbative kx smearing in high-p^ cross sections has the same origin as in 
the 1ow-Qt Drell-Yan cross sections described by the Qr-resummation formalism |I7[| . More re- 
cently, Li [[18] has shown how threshold and transverse momentum resummation may be derived 



from the same parton distribution, defined in transverse momentum space, as in Ref. |§. 

When the conservation of energy is taken into account, however, it is no longer obvious 
whether the inclusion of recoil effects will lead to an enhancement or a suppression, because the 
extra radiation involves a number of competing effects. On the one hand, a substantial kx from 
initial-state radiation allows a softer 2 — > 2 subprocess at the hard scattering, which clearly acts 
toward enhancement. On the other hand, the extra energy of the initial-state radiation drives 
the physical parton distributions to larger x, which may more than make up for enhancements 
in the hard scattering if the distributions are decreasing with x. At the same time, larger x is 
associated with larger threshold enhancements in general. The only way to estimate the influence 
of recoil effects on cross sections is to develop a self-consistent resummation formalism. 

In this paper, we shall take up, and we hope clarify, this general viewpoint. Our reasoning 
is based on a generalization of threshold resummation which, as we have seen, controls singular 
distributions at z — 1. For electroweak annihilation and single-particle inclusive cross sections, 
such contributions are always associated with an underlying 2 — > 2 hard scattering |J. We use the 
2^2 subprocess to define the relevant transverse momentum Qt, whose singularities we resum. 
The recoil we discuss below is always the recoil of a 2 — ► 2 subprocess. Thus, just as for threshold 
resummation, we reorganize a well-defined set of higher order corrections in hard scattering 
functions, always working at leading power in the hard scale, Q, within collinear factorization. 
We do not exclude the possibility of nonperturbative effects, however. Indeed, we will observe that 
nonperturbative corrections arise quite naturally from our resummed expressions. A summary 
of our results, applied to prompt photon cross sections, was described in Ref. [ 19fl . 



Let us offer a few additional comments on nonperturbative effects in these cross sections. 
Nonperturbative effects play a crucial role in the phenomenological description of electroweak 
annihilation cross sections at low Qt, even for a large final state mass scale Q. This is the case, 
even after the resummation of logarithms of Qt / Q that can give a well-defined perturbative pre- 
diction for small Qt- Incorporating nonperturbative effects, of course, requires the introduction 
of new parameters [§, |2(], In each of these cases, the form of nonperturbative corrections 
is suggested by perturbation theory |22|, |23], ^4j]. In contrast, nonperturbative effects (beyond 
fragmentation functions) have not been incorporated in prompt photon and other single-particle 
inclusive cross sections, where there is no need for them at fixed order in perturbation theory. For 
threshold-resummed cross sections the situation is somewhat more subtle, but "minimal" formu- 
lations of threshold resummation allow for a class of purely perturbative predictions |25], |2"6f , with 
no new parameters. Of course, the existence of such a formalism does not by itself preclude the 
importance of nonperturbative effects. In this paper, we develop a perturbative formalism that 
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links both sorts of cross sections, and which is consistent with known results that have suggested 
nonperturbative corrections at measured Qt in electroweak annihilation. Part of our goal is to 
open the door, not only to further perturbative analysis, but also to the study of similarities and 
differences in the roles of nonperturbative corrections in these cases. 

We choose to work in the formalism of collinear factorization because we do not wish to in- 
troduce a new set of phenomenological parton distributions, depending on transverse as well as 
longitudinal degrees of freedom, except where absolutely necessary []. The resulting combination 
of threshold and transverse momentum resummations is at least as technically challenging as 
NLO factorization, let alone QT-resummation, and the new formalism will require some time 
to understand and develop as a practical tool. We therefore do not attempt to draw imme- 
diate phenomenological conclusions in this paper. Instead, we shall concentrate on the formal 
development, and (especially in Appendix A) the theoretical underpinning of these ideas. We 
have attempted to be as explicit as possible in our arguments and in specifying the functions 
whose momentum- dependence controls the set of higher-order corrections that we study. This 
has resulted in a paper of substantial, although we hope not excessive, length. 

We begin in Sec. 2 with a treatment of electroweak annihilation processes, such as Drell- 
Yan and Higgs production, whose singular behavior at vanishing transverse momentum has been 
studied intensively over the years M, M, and which is in many ways the archetype for resumma- 
tion. We show how to introduce threshold resummation consistently at measured Qt for these 
processes. Our approach to resummation is through a "refactorization" of partonic cross sections 
near threshold fTl ^7|. In this discussion, we shall review the refactorizations at the basis of Qt 
and threshold resummations, and define a set of new functions which control singular behavior 
in 1 — z and Qt jointly. These will serve as building blocks both for electroweak annihilation 
cross sections, in Sec. 3, and for single-particle inclusive processes, in Sec. 4. 

Resummation at threshold and in transverse momentum is most often formulated in Mellin 
(N) moment space for the former, and impact parameter (b) space for the latter. Resummed 
logarithms of these parameters exponentiate in the relevant limits, so that the resummed cross 
sections are inverse transforms. In Sec. 3 we resum logarithms of b and N in the electroweak 
annihilation cross sections. We begin by deriving a relation for the hadronic da/dQ 2 dQ T in terms 
of parton distributions, eikonal cross sections for partons, and universal anomalous dimensions. 
We observe that this jointly-resummed cross section determines the pattern of power corrections 
in Q and b that are implied by the behavior of the strong coupling in perturbation theory. 
In particular, we find that in QCD such power corrections appear only at even powers of the 
invariant mass Q and impact parameter b. 

Sec. 4 deals first with prompt photon production, and then with general high-p^ single-hadron 
or photon inclusive cross sections. For the former, we derive the joint resummation applied in 
Ref. flnj, and for the latter we discuss the additional resummation associated with fragmentation. 

Explicit NLL expressions for jointly-resummed exponents in electroweak annihilation and 
prompt photon production are given in Sec. 5, along with a few comments on the source of 
enhancement at NLL. Following our conclusions, we include two appendices. The first gives the 
necessary arguments for factorization and refactorization, and the second gives explicit one-loop 
results for some of the functions that play an important role in the refactorizations of hard- 
scattering cross sections. 



1 This may well be the case for vector boson production at low Qt |20|. Nevertheless, we feel that it is important 
to explore fully the simpler formalism. 
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2. Refactorization for Electroweak Annihilation 



As above, Q denotes the mass of an electroweak final state, such as a vector boson, a Drell-Yan 
pair or a Higgs boson. The cross section da/dQ 2 dQ' T , at measured Q 2 and Qt, given in factorized 
form in Eq. (P, is singular at Qt = 0, order-by-order in perturbation theory. There are a number 
of phenomenological applications of QT-resummation for these singularities |?0 , . We know of 



no simultaneous application of threshold resummation, however, although the cross sections are 
singular as well at partonic threshold, s — Q. 

As pointed out above, distributions that are singular at threshold are smoothed in the physical 
cross section by integration with the parton densities. Nevertheless, there is a good deal to be 
learned by resumming singular behavior from the limit s — > Q 2 , even at fixed, measured Qt 



28fl . To derive a cross section resummed both in threshold and Qt variables, we study partonic 
cross sections a + b — > V + X near threshold, where a and b are partons, and where V denotes 
the heavy electroweak final state V = W, Z, H, etc. We begin by formulating the problem in a 
standard form, through collinear factorization. 

In this section, which is rather technical in parts, we lay the groundwork for our derivation 
of jointly resummed cross sections. We have chosen to present our new formalism in the context 
of a review of existing resummations, for the purpose of motivation, and also to bring together 
a set of results and methods that are somewhat scattered in the literature. It may be helpful, 
therefore, to outline the contents and aims of the subsections that follow. 

We begin (subsection 2.1) by relating the hard-scattering functions that we will resum to par- 
tonic cross sections. In subsection 2.2, we review existing refactorizations for partonic cross sec- 
tions, which have been used to derive resummations for electroweak annihilation at low transverse 
momentum H] and at partonic threshold & Q . We then go on to present a novel refactorization 
that combines the two (Eq. @), and observe how refactorization provides a natural formulation 
of the effects of recoil. The all-orders justifications for all of these refactorizations are presented 
in Appendix A. Eq. (§) involves new perturbative functions, denoted TZi/j. The field-theoretic 
content of these functions is the subject of subsection 2.3, which begins with a review of the 



analogous definitions for light cone parton distributions [^9| , as well as the fixed-energy distribu- 
tions introduced in Ref. ||. Each of the refactorization theorems in subsection 2.2 also includes a 
function that describes coherent radiation, which summarizes the interference between emission 
by incoming and outgoing hard partons. An analysis of coherent radiation is especially important 
for processes in which colored particles emerge from the hard scattering JF], the simplest of which 
is prompt photon production. This interference may be treated in eikonal approximation. The 
analysis of the eikonal approximation to soft gluon radiation is the subject of subsection 2.4, in 
which various eikonal analogs of the densities in 2.3 are introduced. Finally (subsection 2.5), we 
review the use of Mellin and Fourier transforms to isolate hard-scattering functions. The new 
results derived in this section are applied in Sec. 3 to electroweak annihilation, and in Sec. 4 to 
single-particle inclusive cross sections. 

2.1. The hard-scattering function 

Although the hard-scattering function a a b^v m Eq. ([[]) is singular at Qt = 0, these singularities 
may be determined at the same time as threshold singularities at Q 2 /s = 1. To be specific, we 
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shall derive an expression for a a b->v in terms of its moments with respect to z = Q 2 /s, 

oab^v (N, Q 2 t /Q\ Q 2 /fi 2 , fi, a a (ji)) = J Q dz z"- 1 a ab ^ v (z, Q 2 T /Q 2 , Q 2 /fi 2 , //, a a (//)) , (2) 

where for economy of notation, we denote the transform of <7 a &_>y with respect to z {Qt) by its 
argument N (b). The hat refers to its role as a hard-scattering function, from which collinear 
divergences are subtracted. The inverse of the Mellin moment @ is, as usual, 

/ \ rC+ioo (IN / \ 

&ab^v (z,QT/Q 2 ,Q 2 /li 2 ,fM,a s (ijL)) = — z~ a a b->v (N,Q^/Q 2 ,Q 2 /fi 2 ,ij,,a s (fi)) . (3) 

We are able to construct o a b-,v{N) because it also emerges from the factorization of partonic 
cross sections with respect to t = Q 2 /S. Up to corrections due to parton mixing, which we may 
neglect to leading power in the moment variable, we have, from Eq. (fl|), 

<W (N, Qt/Q 2 , A*, O.(A0) = T 777-: U 777-: , /' & . (4) 

v ' (f) a/a (N (j) b/b (N + Jo dQ'dQ 2 T 

The moments, <fti/i(N+ 1, /i), of the parton-in-parton distributions cancel collinear singularities in 
the moments of the partonic cross section, and the right-hand side of this expression is infrared 
safe, order-by-order in perturbation theory when Qt 7^ 0. Our goal now is to determine the 
singular structure of a a b^v at both z — 1 and at Qt = 0. To control these singularities, 
we follow Refs. |5|, |30], [H], |32], 34], and refactorize the (collinear-regularized) partonic cross 
section do a b~*v I dQ 7 'dQ T in this limit. The discussion of the following subsection applies entirely 
to these, purely partonic, cross sections. 



2.2. Refactorization and recoil in the partonic cross section 

To motivate the refactorization appropriate to joint resummation in Qt and 1 — z, it may be 
useful to review the relevant features of the separate resummation formalisms for transverse 
momentum and threshold. We will continue to work in the context of perturbation theory as in 
Eq. because our aim is always to analyze higher orders in partonic hard-scattering functions. 
Each of the refactorizations given below involves the introduction of new parton distributions, 
variously at measured transverse momentum and/or energy fraction. The new functions are not 
to be interpreted as physically-accessible distributions. Rather, they are perturbative constructs 
useful for the analysis of the hard-scattering functions of Eq. (Hf). 



In the formalism of Ref. p5fl , the measured-Qr cross section is written as a convolution of 
(parton-in-parton) distributions T > i/j(x,'k), at fixed parton transverse momentum k, and light- 
cone momentum fraction x, along with an additional, eikonal function £/ C( j(q) that describes 
coherent soft-gluon emission at fixed transverse momentum, 



da. 



ab^V 



dQ 2 d 2 Q T 



a fd^v(Q 2 ) h< cd\ a s(Q)) J dx a d 2 k a V c/a (x a ,k a ,Q) J dx b d 2 k b V d/b {xbM,Q) 

cd 

x J d 2 q £/ cd (q//i, a a (n)) 5{Q 2 - x a x b S) 5 2 (Q T + k a + k b + q) + Y kt , (5) 



where a^ is the Born total cross section for the process, for example, cr^ * = 4na 2 /3Nc Q 2 , 
with Nq the number of colors. The remainder, Y kt , does not diverge as a power at Qt = 0. Note 
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that because the functions V%/j are defined at measured k, no factorization scale is necessary, 
although the distributions still depend on the overall momentum scale Q. The variable /i in the 
arguments of U is therefore a renormalization scale. The additive convolution in this expression 
implies that the cross section breaks up into a product under a Fourier transform to impact 
parameter (b) space ||. The function h^\a s (Q)) = l + 0(a s (Q)) absorbs hard-gluon corrections 
that appear in coefficients of 5 2 (Q,t)- The combination ftA kt <* is a truly short-distance function, 
dominated by lines off-shell by 0(Q 2 ) | 36fl . In contrast, the full hard-scattering function in Eq. (|]) 



in general contains lines that are off-shell only by 0(Qj,). This hierarchy of perturbative scales 
is characteristic of resummation. We shall use the term "short-distance" to refer specifically to 
functions that depend only on the largest scale, Q in this case. 

The refactorized cross section for threshold resummation, with integrated Qt, has many 
of the same features. Now, however, the parton-in-parton distributions ipi/j(x) are defined at 
measured fraction of the energy of parton j (in the center-of-mass frame for the hard scattering), 
rather than light-cone fraction, as is the new eikonal function U c d(w s ), with total energy w s Q 
for soft gluon radiation into the final state. Working to leading power in 1 - Q 2 /S leads to 
important simplifications. First, the nondiagonal parton-in-parton distributions, ip c / a (x, Q) begin 
at order a s with the emission of a soft fermion (not a pair) into the final state, which results in a 
suppression of order 1 — x in the distribution, and of 1 — Q 2 /S in the cross section || . To leading 
power, therefore, we may neglect parton mixing, just as at leading power in the moments, Eq. 
(H). The refactorized expression is f| 

d0 ^Q2 = ^ a ab^v(Q 2 ) h ab\ a s(Q)) J dx a 1p a/a (x a , Q) J dx b 1pb/b(x b , Q) 

x J dw s U ab (w s Q/fx, a s (n)) 5 (l - Q 2 / S - (1 - x a ) - (1 - x b ) - w s ) + Y th , (6) 

where Yth is nonleading by a power of 1 — Q 2 /S. Even though transverse momenta in ip are 
integrated, the phase space for radiation is finite for fixed parton energy, and \i again denotes 
the renormalization scale. An explicit definition of ip c / a as a matrix element will be given below. 
The remainder, Yth, does not diverge as a power of 1 — Q 2 /S at threshold. The short-distance 
function = 1 + 0(a s ) organizes infrared safe coefficients of 5(1 — z) in this case. 

It is most natural to analyze the cross section near threshold, Eq. ([]), in terms of a Laplace 
transform, / dr exp[— iV(l — r)], with r = Q 2 /S. For iV large, we can readily relate this Laplace 
transform to the Mellin moments in Eq. Q). This follows because generally, 

e -iV(i-o^ e iV ) (7) 

with corrections that are suppressed by a power of N, and because in Eq. (|6|), 

1 - (1 - x a ) - (1 - x b ) -w s ~ x a x b {\ -w s ) + ([l - Q 2 /Sf 

The Laplace moments of Eq. (|6]) are therefore equivalent to its Mellin moments to leading power 
in N, and hence in 1 — Q 2 /S. 

The close correspondence between the factorizations at low Qt and near threshold makes 
it rather natural to combine the two. We therefore propose a convolution at fixed transverse 
momentum and energy fraction: 

dQ^Q T = ^ a ^v(Q 2 ) h ab( a s(Q)) J dx a d 2 \t a H a / a (x a ,\L a ,Q) J dx b d 2 k b K b/b (x b ,k b ,Q) 
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a 




Figure 1: Leading region for electroweak annihilation in cut diagram form. The vertical line 
represents the final state that includes an electroweak boson, labelled Q. The subdiagrams J a , 
Jb, U a b and H incorporate, respectively, on-shell lines with momenta collinear to p a , lines with 
on-shell momenta parallel to pb, lines with soft momenta, and lines off-shell by order Q. 

x J dw s d 2 k s U ab (w s ,Q,'k s ) 5(1 -Q 2 /S - (1 - x a ) - (1 - x b ) -w s ) 

x 5 2 (Q r + k a + k b + k s ) + Yi . (9) 

The short- distance function h^(a s ) is again an infrared-safe series in the running coupling, 
which begins with unity at zeroth order, and which absorbs, in this case, the coefficients of 
5(1 — z)5 2 (Qt) at one loop and beyond. The remainder Yj is free of power singularities at 
Qr = at leading power in 1 — Q 2 /S. As in threshold resummation, only flavor-diagonal hard 
scatterings contribute at 0[1/(1 — Q 2 /S)}. It is important to note that in terms that are not 
singular in Qt, this leading power emerges only after integration over Qt- This is because at 
fixed energy (1 — z)Q, the phase space in Q T behaves as: J Q dQ\ = (1 — z) 2 Q 2 . 

Eq. (|S]), and indeed each of the refactorizations discussed above, may be represented as in 
Fig. 0. In the terminology of Ref. || and Appendix A below, Figure [I] represents the general 
"leading regions" in momentum space for this cross section. The subdiagrams J a ^ include lines 
collinear to the incoming partons, H lines off-shell by order Q, and U soft radiation. 

The refactorizations of Eqs. fl5|) and @ themselves define the concept of recoil that we will use 
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in this paper. The short- distance function is computed with on-shell external momenta, 

collinear to the incoming lines. All unintegrated transverse momentum dependence is contained 
in the generalized parton densities V in (|5|) and 1Z in Eq. ([5]). The dependence of highly off- 
shell lines on the transverse momenta k a and kf, of initial-state partons is to be absorbed into 
higher orders of the short-distance function, by the usual methods of collinear factorization. On 
the other hand, in both transverse momentum and joint resummation, we retain the kinematic 
linkage of the partonic transverse momentum with the electroweak final state. This is what we 
shall mean by including recoil effects. 



2.3. Matrix elements 

The refactorization theorems above, and the resummations derived from them, involve a number 
of new functions. We now give explicit definitions for the various parton distributions, V, TZ 
and t/>, when the incoming partons are quarks, as well as for the eikonal functions U . Gluonic 
distributions can be defined similarly, following Ref. 



The parton densities lZ a / a and the eikonal functions Ufj, defined at fixed energy and transverse 
momentum are, like Eq. (|9]) itself, straightforward variations of functions identified for the Qt 
and threshold resummed cases. The prototype for these expressions is the partonic light-cone 
distribution, written as 



ji,e) = -L /^e"^ + </(p)|ff / (Au) 7 -tig/(0)|/(p)> > (10) 

where \x is the scale at which the product of quark fields, which are connected by a lightlike 
separation, Aw M , u 2 = 0, is renormalized. An average over colors and spins is included in the 
definition. In this expression, we have suppressed an ordered exponential, $^ (A, 0;0), which 
we shall also refer to as a nonabelian phase line, of the gauge field along the light cone vector 
between the quark fields, in the notation 

^(X 2 ,Xi',X) = Pexp(-ig J*' drj + X) \ . (11) 

Here the gauge field is a matrix in the representation of parton /. In momentum space, these 
operators correspond to eikonal lines. Equivalently, we may define the matrix element ( |10"D 
in u ■ A = gauge. The perturbative MS distribution, computed as a power series in a s (//), 
is independent of the momentum it is a "pure counterterm" , that is, a series of poles in 
e = 2 — D/2 in D dimensions, a dependence which we exhibit among its arguments. 

We may regard perturbative MS distributions as defined by their evolution equations, which 
in moment space are 

y 2 -^ h/f(N, n, e) = 7// (N, aM) 4> ff (N, //, e) + 0(1/ jV) , (12) 

with 7jj (N, a s , e) = Jq dxx N ^ 1 Pff(x, a s ) the moments of the splitting function for flavor /. As 
usual, up to corrections of order 1/N, we may neglect flavor mixing. A very useful explicit 
form for the MS distributions is found by solving this equation, with the boundary condition 
4>(N,n = 0,e) = 1, 



4>f/f( N ^^ e ) = ex P 



2 ^,,'2 



djj: 



1 ff (N,a 8 {fj?)) 

h 1 



(13) 
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This expression is meaningful for the collinear-regularized distribution, defined for D > 4, or 
equivalently e < 0, because of the e-dependence of the strong coupling. The p! integral is 
regulated by reexpressing a s (p') in terms of the strong coupling evaluated at the fixed scale \i: 
a s (/i) = (p' '/ p) 2e a s (p') + . . .. The MS prescription then consists of reinterpreting the upper limit 
as fMp ln(47reg 7 ), with 7^ the Euler constant. We shall not generally exhibit this modification 
below, nor indicate explicitly the e-dependence of the coupling. 

To leading power in N, jff is found from moments of the expansion 



P ff {z) = A f (a s ) 



1 



ll-z 



B f (a s )S(l-z) + 0([l-z] ) . (14) 



The plus distribution appears only as a power series in a s times [1/(1 — z)]+. It is worth noting 
here that the form of Eq. flUD, with no explicit powers of ln(l — z) in the plus distribution, 
is required by collinear factorization, and is not an additional assumption [|37|]. We shall see 



this result emerge below in Section 3. A similar observation was made recently by Albino and 
Ball [0. 



For NLL expansions, we will need the anomalous dimensions A a (a s ) to two loops. For flavor 
a, they are given by the familiar expansion, A a (a s ) = Y, n ( a s/^) n A^\ withP] 



4 1} = c a 

4 2) = \c a K = \c a 




(15) 



where C q = Cp, C g = Ca- To lowest order, which is the accuracy necessary for NLL, we have 
B a = (a s /ir)B^\ where Bjp and Bjp are given by 

BV = \c Fi BV = S±, (16) 

with (3q = HNc/3 — 2rif/3, the lowest-order coefficient of the QCD beta function. 
Matrix element representations of the functions Vf/f are similar |29 , 



1 r rl\ rPh 

V f/f (x,k,p.n,e) = —J — ^e-^^ b - k (/( P )|g7(0 + ,A,b) 7 - Mg/ (0)|/( P )). (17) 

This matrix element is defined in an axial, n ■ A = gauge, which is how it acquires p ■ n 
dependence. It also requires collinear regularization in perturbation theory. 

The densities ipf/f are the distributions of quarks of fixed energy p = Q/2x ~ Q/2, in the 
center-of-mass of the produced pair, while lZf/f(x, k, p°) are distributions in energy and transverse 
momentum k. The external line is an on-shell quark of four-momentum p M = (Q /x^/2)v fM , with 
v 2 = 0. The inclusive energy distribution ipf/f is then given by 

Po f d\ _ iXxpo 



^ f/f (x,2p ,e) = J ^ e-^° (/(p)|ff/(An) 7 • « ?/(0)|/(p)> , (18) 

where qj{x) is the field for flavor /, is the light-cone unit vector opposite to u ■ v = 1, 
and is the unit vector in the time direction, = (1,0). Following Ref. 0, we evaluate the 

2 The function A a is proportional to r cusp in Refs. @ 0. 
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matrix element fljjp in n ■ A = A = gauge in the center-of-mass frame, which turns out to 
be convenient for calculational purposes. Were ipf/f defined in a spacelike axial gauge, it would 
differ only by finite corrections. The operator product separated by a timelike distance requires 
no new renormalization. Correspondingly, the functions 7Z q / q may be defined as matrix elements 
by 

n f/f (x,k,2p ,e) = ^^/|) ^ e"^ b ' k (miqfiXn + bh.uqfmm), (19) 

again evaluated in A = gauge. In these expressions, and in the remainder of this section, we 
suppress dependence on the renormalization scale, which we take equal to Q = 2po. We will 
return to the choice of renormalization scale later. 



2.4. Eikonal functions and factorizations 

Near partonic threshold, all radiation is soft, compared to the hard scattering function. It is thus 
natural to study the eikonal approximation for the cross section and for the factorizations that 
characterize the dynamics. The discussion below follows Refs. 

The eikonal cross section is built from ordered exponentials, i^j (0, — oo; X), of the gauge field 
in the group representation of the incoming partons, extending from minus infinity to the point 
of annihilation, in the notation of flTTp. We introduce a product that represents the annihilating 
combination of two nonabelian phase operators: 

W^(X) = 4f ) (0,-oo;X) 4 c) (0,-oo;X) , (20) 

where for quarks, c and d may carry different flavors, as in the case of u + d — > W + . From the 
operators W*"^, we define an eikonal cross section at fixed energy and transverse momentum, 
which represents the QCD radiation generated by the annihilation of the two incoming color 
sources, neglecting recoil, 

dX d 2 b 



fx (eik) ^ O k e) - Of— e -iA?Q/2+ib.k 



2tt (2tt) 

x — Tr (01 T [W (cj) (0) f ] T \w {cd) (Xn + b)] 10) . (21) 
d(c) Lin 

As above n M = 5^q, so that Xn + b represents the vector with time component A and transverse 
components b. In the matrix element, T represents time order and T anti-time order. The trace 
is over color indices in the representation of parton c. d(c) is the dimension of this represent at ionP]. 
Because the velocities (3 and f3' of the incoming lines are lightlike, this cross section has collinear 
singularities, and must be regulated. Infrared divergences, however, cancel in the sum over final 
states @]. 

To organize collinear singularities, we introduce eikonal parton distributions, which approxi- 
mate the radiation at fixed energy and transverse momentum from an energetic, lightlike parton, 

x Tr (0|4 a)t (0,-oo;0) ® { "\o, -oo; y n + b)|0) , (22) 



3 The eikonal cross section denned here is normalized to S(£,)S 2 (k) at zeroth order. The average over the colors 
of the physical incoming partons will be absorbed into a separate overall factor. 



11 



computed in the n ■ A = gauge, just as 72-///, Eq. (|19D . Similarly, by analogy to Eq. (|10D , we 
can construct an eikonal distribution at fixed light-cone momentum fraction, 

9. [ dX -i£QVV2 



(eik) 

y/2 d(a) J 2tt 

x Tr (0|$^ a)t (0,-oo;0) $£°(0, -oo; Au)|0) , (23) 

which as usual requires renormalization of its ultraviolet divergences, and regularization for its 
collinear divergences. As in Eq. (|To|) , we omit the ordered exponential in the opposite- moving 
light cone direction, between and Xv. Also like the MS distribution, 0///, 0i° lk - ) is a pure 
counterterm, and is independent of the momentum scale Q and of the direction of (3. It is also 
flavor-independent among quarks and antiquarks, differing, of course, for gluons. Note that £ in 
0^ eik ) plays the role of 1 — x in 0///. That is, we fix the light-cone component of the emitted 
radiation, since the eikonal line does not have a definite initial-state momentum. 

Purely virtual diagrams in both 7e (eik) and <M eik) enter as overall factors, which can be used 
to normalize these functions. We choose to define the virtual contributions by the requirements 
that 

f l dw [ d 2 k Kf k \w, Q, k, e) = /' d£ 0i eik )(£, /i, e) = 1 . (24) 
Jo J Jo 

These conditions ensure that both functions are sums of plus distributions in terms of the variables 

( or w, integrated over the interval from zero to unity. This choice does not affect the N- 

dependence of the functions at all, but ensures that factorization does not introduce spurious 

collinear singularities. This condition also enables us to define an evolution equation for the 

eikonal light-cone distribution of the form of Eq. fll2|) , with solution 

d^' 2 (eik) 



4>f k) (N,fi,e) = exp 



r 1 



(25) 



which differs from fll~3|) only in the eikonal approximation to the anomalous dimension. The 
eikonal anomalous dimensions, 7// k ' ) , are found from the plus distributions of the splitting func- 
tions, when written as in Eq. (|T4|), subject to the normalization condition (|24|) . To leading power 
in N, the moments of the eikonal distribution in D dimensions are given by 

dfi' 2 



(f)f lk) (N,fi,e) = exp 



Jo 



(26) 



where we define 



N = Ne yE . (27) 

The discussion on the dimensional regularization and MS definition of <f> a / a {N, fj,,e) given after 
Eq. (0) applies as well to its eikonal analog in Eq. (p5[). 

Essentially the same arguments (see Appendix A) for the joint factorization of the partonic 
cross section are valid for the eikonal cross section, cr ^- We may therefore factorize cr^| lk ' ) , in 
terms of energy and transverse momentum distributions, as in Eq. @, 

^ lk) (£,Q,k,e) = 

J dw c d 2 k c 7^ cik > (w c , Q, k c , e) J dw d d 2 k d TZ d cik) (w d , Q, k d , e) 

x J dw s d 2 k s U cd {w s , Q, k 8 ) 5(£ - w c - w d - w s ) S 2 (k + k c + k d + k s ) . 

(28) 
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Alternately, we may factorize the eikonal cross section in terms of eikonal light-cone distributions, 



cS^Q^e) = / dw c (f)^ lk \w c ,fi,e) J dw d <pf k) (w d , fi, e) 

x J 'd? 5(Z-w c -w d - O *S k) (r , Q, k) , (29) 

where partons c and d are implicitly in a color singlet state. This is the eikonal approximation to 
Eq. (HD for incoming partons c and d, identifying Wi ~ 1 — x iy £ ~ 1 — Q 2 /^ and £' ~ 1 — Q 2 /s. 
In Eq. (p8|) and (|29|), respectively, the 7?.'s and 0's absorb collinear singularities associated with 
soft gluons. The remaining functions U c d and are then infrared safe. 

All of the refactorizations in Eqs. (H), @ and (||) involve functions that are gauge-dependent. 
We have already noted the gauge- dependence of the functions 1Z and TZ^ eik ^ above. The soft 
function U ab also inherits gauge dependence through TZ^ clk \ To the extent that the factorization 
formulas are valid, however, all gauge dependence is guaranteed to cancel in the cross sections. 
In Appendix A we study the theoretical basis of these refactorizations; for the purposes of the 
following discussion, we accept their validity. 

2.5. Transforms and the soft function 

As usual, the refactorized cross sections are displayed most conveniently in terms of their ap- 
propriate transforms, in this case, Laplace and Fourier. Transforms that we will need below 
are: 

H f/f (N,bQ,e) = J™ dxe- N{1 ~ x) J d 2 ke- lh ' k n f/f {x,Q,k,e) 

_ poo p 

7^ eik) (N, bQ, e) = J dw e~ Nw J d 2 k e~ ibk TZ { a c[k) (w, Q, k, e) 

U cS (N,bQ) = J dwe~ Nw J d 2 ke~ ihk U cd (w s ,Q,k) 

af\N, bQ, e) = jT d£ e"^ / d 2 ke~^ ag k) (£, Q, k, e) . (30) 

These transforms simplify the double convolutions of the partonic cross section, Eq. (|9[). On the 
other hand, the moments of the hard scattering functions of Eq. ([[]) are determined from the 
partonic cross sections via Eq. (|). In this way, we find 

ta^v(N,Q 2 T /Q 2 ,Q 2 /n 2 ,n,a s (ri) = a%l v (Q 2 ) - — L 

x J ^ e* bQ - TZ a/a (N, bQ, e) K b/b {N, bQ, e) U ab (N, bQ) , (31) 

where we define a combination of short-distance function and Born cross section as 

^ly(Q 2 ) = tt ^l v {Q 2 ) h^Mn)) . (32) 

The factor 7r relates the azimuthally symmetric da / dQ 2 d 2 Q,T to da j dQ 2 dQ\, while the Born 
cross section af b \ v absorbs the color average for the initial-state partons, referred to above. In 
Eq. ([y]), we have approximated N + 1 by iV in the arguments of the light-cone distributions, as 
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is acceptable to leading power in N. Again, we suppress the renormalization scale fi = Q; the 
explicit \i here has the interpretation of a factorization scale. 
Another useful form of Eq. ([[!]) is 



Oil 



oab^v (N, Q 2 T /Q 2 ,Q 2 /fi 2 ,fi, 

= *%lv(Q 2 ) J 7S2 eib ' QT ^/a(N, b, Q, fi) c b/b (N, b, Q, /,) , (33) 
where the functions c and their eikonal analogs are defined by 

- (ATUn , n f/f (N,bQ,e) \O ff (N : bQ) r ' 



1/2 



^ eik) (iV, bQ,e) [%(i V,bQ) 



gpWb.g,/*) = — ' ^ — ^— . (34) 



The c's will appear as building blocks in direct photon and other cross sections with color flow 
into the final state as part of the hard scattering. 

The eikonal cross section has many of the same properties as its partonic counterpart. Mo- 
ments of the eikonal hard-scattering function at fixed k are found from (^9|): 



a (eik) (N O k u) ~ &(X,Q,k,e) 

& ; (iV,/x,e)^ '(N,fj,,e) 



where to avoid unnecessary clutter in our notation, we identify the transforms of the functions a 
only through its arguments. Then, using the eikonal transforms in Eq. fl3"0| ) in the eikonal joint 
convolution fl28|) , we have 



0( cik) (iV,/i,e) 0f ik) (iV,/i,e) 

d 2 h 

(2tt)' 



x / e* bk n£ ik \N, bQ, e) ^ cik) (iV, bQ, e) U cd (N, bQ) , (36) 



with the same coherent function [7 c j. 

For completeness, and for reference below, we observe that the eikonal function Uij at mea- 
sured energy and transverse momenta may be defined by its transforms, through 



*% k \N,bQ,e) 
Hr K) (N, bQ,e) Tl { f k \N,bQ,e) 



UUN, bQ) = ; , (37) 



where collinear singularities cancel in the ratio. The other forms of the soft function in Eqs. @ 
and (|5|) differ only in the components of the total final-state momentum that are fixed. Note 
that this expression is essentially a rewriting of the refactorization for the eikonal cross section, 
Eq. ©. 

The behavior of the hard scattering function at large moment N and impact parameter b 
may be studied either in terms of U and the distributions 71, or, as we see in the next section, 
by relating the partonic and eikonal functions given in Eqs. fl3~T|) and (|36|). In the remainder of 
the paper, we apply this formalism to derive our jointly resummed cross sections. 
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3. Joint Resummation for Electroweak Annihilation 



As pointed out in Refs. |27|, [43|, |44}|, color-singlet cross sections with symmetric phase space expo- 
nentiate at high moments that force the phase space to an "elastic" limit, where only soft gluon 
radiation is allowed. This is the case for doubly-transformed cross sections, and for the singly- 
transformed cross section in threshold resummation |5], Q . The elastic limit is naturally associated 
with the eikonal approximation. As we now show, the full leading-power N- and 6-dependence 
of electroweak annihilation cross sections can be deduced for quark- ant iquark annihilation and 
gluon fusion directly from eikonal cross sections. 

3.1. Partonic and eikonal cross sections 

Near threshold, that is, to leading power in N, all real-gluon emission in the partonic hard- 
scattering function, Eq. ([H]) may be treated in eikonal approximation, Eq. (0). Since the 
function U is the same in the partonic cross section and its eikonal approximation, the difference, 
for fixed b, resides entirely in the parton distributions, and we have 

K a/a (N,bQ,e) n b/b (N,bQ,e) = n^(N,bQ,e) 7^ eik) (iV, bQ, e) 

4> a/a (N,}x,e)^ /b (N,Li,e) ~ 4> { : ik \N, t x,e)4> < t k) (N, t ,,e) ' 

where V ab (Q,fx) is an overall factor, entirely from virtual corrections, which are of the same 
graphical form in 1Z and <fi. The function V a b{Q, aO is therefore independent of b and N to order 
N°; its Q and /i dependence may be determined as follows. 

As shown explicitly in Eq. (|T3| ) and (]25|) above, light cone distributions and their eikonal 
approximations in the MS scheme are fully determined by the splitting functions. In moment 
space, the leading power in N comes entirely from the transforms of [1/(1 — z)] + and 5(1 — z) 
contributions to the splitting functions of Eq. ([14]). The former, which can only arise from the 
combination of real-gluon and virtual corrections, are fully represented in the eikonal distribu- 
tions, (fi( mk \ As a consequence, to leading power, the ratios (f)^ lk \N, fx) /(p a / a (N, fx) only receive 
contributions from the left-over 5(1 — z) terms in the splitting functions, 

~l — ~r T rr (AT r=exp - —^ w [B a (a s (fi)) + B b (a s (fi))]\, (39) 

with fx p the factorization scale, with explicit dimensional regularization, and with the functions 
B a given by Eq. (|I"6|). The collinear divergences in this expression cancel in the ratio in Eq. (3S). 



The ratio of the ^-functions must thus take the same form at leading power in N, but with 
an upper limit on the //-integral given by the renormalization scale in TZ, which, as above, we 
choose to be Q, 

exp - / [B a (a s (n')) + B b (a s (ix'))} . (40) 



K a/a (N,bQ,e)n b/b (N,bQ,e) r 1 Jo ix' 2 +^ 

As a result, the uniquely determined form of V ab (Q,Lx) is 



f rQ 2 du' 2 1 
V ab (Q,n F ) = exp | [B a (a s (fx')) + B b (a s (fx'))] \ . (41) 
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For electroweak annihilation, of course, B a = B b , but the result is quite general. 

We are now ready to combine Eqs. (|3lD for the partonic hard-scattering function, Eq. (^) for 
its eikonal approximation, and Eq. (f4T| ) for the ratio V, to derive an expression for the refactorized 
electroweak cross section that we will study below. The result is: 



(N,Q 2 T /Q 2 ,Q 2 /fj, 2 F ,fi F ,a 8 (fi) 



^iv(Q 2 ) exp [B a (a s (fi')) + B b (a s (fi'))] 




x 



(,2h e*><*&W(N,b,Q t/ i). (42) 



(27T) 



Equivalently, the hadronic cross section, given as an inverse transform from moment space, takes 
the form 

d °^ L ~ E ^lv(Q 2 ) exp {/f^ [B a (aM) + B b (a s ^))]\ 



X Ic^i ^ a/A ( N + lj ^f)^/b{N + 1, y. F ) t 



c 2ni 

X /(2^ e ' b ' QT "S k)(iV ' b ' g ' /i) ' (43) 

with fip the factorization scale and r = Q 2 /S. Corrections are order 1/N in the hard scattering. 
The sums over a and b include quarks, antiquarks and gluons. Thus, for the case of electroweak 
annihilation, a direct examination of the eikonal cross section will determine the large-moment 
and impact parameter behavior of the partonic hard-scattering functions. Here and below, C is 
a contour to the right of iV-plane singularities in the various transform functions. 

Eq. (|43|) reduces the computation of the cross section in MS scheme to the determination of 
the eikonal cross sections. Alternately, one may treat the factorized functions ip and U separately, 
applying renormalization-group arguments J5|, [J . The result must be the same. In the remainder 
of this section, however, we shall examine the eikonal cross section directly. 



3.2. Exponentiation 

We employ a number of general features of eikonal cross sections, where the graphical analysis 



introduced in Refs. |]43 , |4q , |4q| is particularly helpful. Moments of cross sections in the eikonal 



approximation exponentiate at the level of integrands, with exponents given by the moments of 



a set of graphical functions that Gatheral termed "webs" . Webs can be generated uniquely 
from cut diagrams in eikonal cross sections order-by-order. They are defined both in terms of 
graphical topology (irreducibility under cuts of the eikonal lines) and color structure. The lowest- 
order web is simply a single gluon exchanged between the lines. Beyond lowest order, each web is 
itself a cut diagram, and can be integrated over the momentum, k, that it contributes to the final 
state. A very useful additional feature of webs is that at fixed k they have no overall ultraviolet 
divergences. 

Quite generally, then, the joint moment and impact parameter dependence of the eikonal 
cross section may be expressed as 
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XWab k -W~" ^ asi ^ £ ) ( e_JV(fe ' ft/Q)_ib ' k - } > ( 44 ) 

where u> a ft represents the web at fixed total momentum fc^ 1 , where n M = 5^0, and where for conve- 
nience we choose the factor Q±-2e to be 27r 1_e /r(l — e) ~ 27r(7re 7J3 )~ e , equal to the dimensionally- 
continued transverse angular integral. In this form, the single-gluon emission contribution to the 
web is normalized to be 



lX \k) = Un^e-^Y -L 5 + {k 2 ) . (45) 

7T ^ ' rCrp 



That the overall coefficient is twice the one- loop term in the function A a (a s ) in Eq. ( |l"5"D is not, 
of course, a coincidence, and we shall see below how this relation arises. The e-dependent factor 
matches the MS collinear subtraction. In Eq. ([0]), we explicitly limit the phase space for gluon 
emission by the plus and minus momenta of the annihilating partons at partonic threshold. As 
in the case of the eikonal distributions TZ^ eik ^ and 0( eik ) in Eq. Q2"4]), dependence on the choice of 
upper limits is exponentially suppressed in N for real gluon emission, but must be specified to 
set the scale of virtual corrections. We determine purely virtual corrections in the eikonal cross 
section by demanding that it be normalized to unity at b = 0, iV = 0. This choice ensures that 
the perturbative expansion of o-( eik ) is a sum of plus distributions. 

The fc-dependence of w a b(k, (3, ft') in Eq. ([44]) is strongly constrained by the invariance of the 
web under rescalings of the light-like eikonal velocities, (3^ and (3'^. We then observe, based on 
the lack of overall UV divergences in the webs, that they obey 4 

w ab [k 2 , k j-pT ,t* 2 ,<* B (tJ>),eJ = . (46) 

We know even more about webs, because the purely virtual web is the logarithm of the lightlike 
eikonal form factor, discussed extensively in treatments of the Sudakov form factor in QCD 



40| , f42| . From these investigations, and by comparison to resummations for the Drell-Yan cross 
section [[| |(|, we learn that the webs may have at most a single overall infrared divergence, 
coupled with a single overall collinear divergence. Additional logarithmic singularities can arise 
only through the renormalization of subdiagrams. As a result, at fixed values of kr, relative to 
the axis determined by (3 and (3', the integral of the web over k 2 is finite. The web integrals are 
ultraviolet divergent for kj- — > oo and are collinear and infrared singular at k^ = once k 2 is 
integrated. 



3.3. The exponent 

The exponentiated eikonal cross section contains a considerable amount of information, which 
follows from the properties of the web described above. We use the azimuthal symmetry of the 
web Wat, in Eq. to organize the transverse and light-cone integrals of k into the form 



4 ik) (iV,bQ,e) = exp 



d 2 ^kn 



2c 



Q 







dk 2 w a b (k 



2 k 2 T + k 2 ,fi 2 ,a s ([i), 



4 We note that to derive this relation, the webs should be renormalized appropriately in terms of their external 
eikonal lines. These renormalization factors vanish in Feynman gauge, and we shall ignore them below. The 
overall result for the eikonal cross section, of course, is gauge independent. 
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V 



Q/V2 



dk + - N ^\h + ^v 



(kl +k 2 )/V2Q 2k + 



In 



Q 2 



h 2 -I- lc 2 



(47) 



where the light-cone variables refer to the frames for which (3t = /3' T = 0. For large N, the k + 
integral is well-approximated (up to exponentially-suppressed contributions) by a Bessel function, 
and we find 



a£ k \N,bQ,e) 



exp < 2 



^l-2e JO 



e -ib k T R 2N 



h 2 -A- h 2 

IX rp — p A, 



-In 



dk 2 w ab (k 2 , k% + k 2 , /j 2 , a s (n), e 



Q 2 



h 2 -I- h 2 



(48) 



In view of our comments above about the infrared sensitivity of the web, we are particularly 
interested in the limit that k\ + k 2 vanishes. From the behavior of Kq(x) for small values of its 
argument, 

K Q {x) ln(xe 7B /2), (49) 

we see that the momentum dependence of Kq cancels the logarithm in the collinear limit, leaving 
a factor ln(iVe 7E ). This remainder generates a collinear logarithmic singularity at k\ = 0, which 
is canceled by the moments of eikonal distributions, as we now show. 

To construct the eikonal hard-scattering functions, a^ elk \ given by Eq. (|35|) , we substitute Eq. 
fl48| ) for the eikonal cross section, and the explicit expression (p5|) for the eikonal distributions, 
into the Fourier transform of cr( eik ) to derive 



°ab <{N,b,Q,fi) = exp 



d 2 - 2t k 



T 



2c 



dk 2 w ab (k 2 , k 2 T + k 2 , jj 2 , a s (fi), e) 



ibk 



T K I 2N* 



h 2 -I- h 2 



-In 



Q- 



h 2 -I- h 2 



d=a,b 



(50) 



where we have relabeled the variable \j! in Eq. (|25|) as kr, and where N is defined in Eq. (f27|). 
Eq. ( |50"D has a lot in common with standard resummations in logarithms of N and b, although 
it still includes an extra integral over k 2 . It may be further simplified, however, using properties 
of the webs. 

We continue by rearranging Eq. (|50"D into a form that isolates double logarithmic behavior, 



a% lk) (N,b,Q,n) = exp 



d 



2-2ei 



r<4 2 -k\ . 

J o dk 2 w ab [k 2 ,kT + fc 2 ,/i 2 ,a s (/i),e 

q j--\Qji + wr- lnN E Adic " ikT,) 



-ib k T 



^1-2<e 

2Nk T \ , (k T 
o I I + In I — 



x exp < 2 



d 2 ~ 2£ k T 



e -ibk T I ^ 2J y 



1,2 I h.2 



dk 2 w ab (k 2 , k\ + k 2 , /i 2 , a s (/i) 
2iVfc T 



Q 



In 



p _|_ L2 

IX rp — p /\, 



.(51) 



Let us deal with the two exponentials on the right-hand side of this relation in turn. The 
first exponential in Eq. ( |51~1) begins at the leading logarithm (LL). We have seen that the webs 
contain no internal collinear or infrared divergences. Also, because the web requires no overall 
ultraviolet subtraction for kx fixed, the k 2 integrals of w a b in Eq. (^) converge on a scale set by 
kj<, independent of N, b or Q. At the same time, the factorizability of the eikonal cross section 
requires the cancellation of all singularities at kx = in Eq. (|5l|) . We may thus formally expand 
the integral of the web over k 2 in inverse powers of Q 2 , with a leading coefficient that behaves 
as \jk\ for kp — > 0, which must match the collinear singularity of the subtraction: 



/' 

Jo 



dk 2 



W a b 



2 u2 

rvj-i 



k 2 ,jj 2 ,a s (/j), 



A a (a s (k T ))+A b {a s (k T )) 

(k 2 y-< 



Aab (a s (k T ),k T ,Q) 



(52) 

where the function A a b behaves as (k^/Q 2 ) for kj< — > 0. In this expression we have used the 
renormalization-scale invariance of the webs, Eq. ([46]), to set the scale of the coupling at hp, 
which is the only remaining kinematic variable. Given this result, the kp integral in the first 
exponential of (|5T|) is seen to be finite, and we may remove the dimensional regularization on 
this factor. Leading threshold logarithms in the perturbative expansion of the exponent, for 
example a™ ln n+1 N, are generated by the explicit logarithms of k? and N, which multiply the 
\jk\ behavior isolated in Eq. (|52|). 

In the second exponential on the right-hand side of Eq. (|5T|), the term in square brackets 
behaves smoothly for k\ ~ k 2 — > 0, as well as for k? — > with k 2 fixed. As a result, this 
term has a finite e — > limit, and begins with next-to-leading (NLL) logarithms, for example, 
a™\n n N. With MS eikonal distributions, however, even NLL logarithms are absent, because at 
leading order the web w a \, is proportional to S(k 2 ) (i.e., one-gluon exchange), which vanishes in 
this factor. 

Using Eq. (|52|) in (0), and setting the number of dimensions to four, we derive an expression 
for the resummed cross section in transform space, 



a a f\N,b,Q^)=ew D a T'(N,b) 



.(eik). 



exp 



(53) 



where the leading iV and 6-dependence (LL and NLL) is entirely contained in the exponent 
E^\N,b,Q,!A F ) 
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(54) 



The second term on the right accounts for the difference between the physical scale Q and the 
factorization scale [ip. 

The factor exp[D a b] in Eq. (^) contains corrections in the form of an infrared safe expansion 
in a s plus NNLL and nonleading powers in iV and b, 
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(55) 



These, rather elaborate, expressions are accurate to all logarithms in N and b, and implicitly 
contain as well the structure of power corrections [^, P3| , which we hope to study in future work. 
Here, we only note that the expansion of the function Kq(x) for small x contains, up to a single 
logarithm, only even powers in x. This simple observation is enough to ensure that the threshold 
resummation of perturbation theory to any order implies the presence of even powers of Q^ 1 
(and b) only |22], |24fl . In this paper, we shall regard the above results as the starting point for 
a joint NLL resummation in N and b for electroweak annihilation, and for high-p^ photon and 
hadron cross sections. As noted above, the entire NLL result is associated with E^, although 
D a b may contribute beginning at NNLL. 



3.4. Hadronic cross sections 



The explicit form of the jointly resummed cross section is now found by inserting Eq. (|53| ) in 
(|43|). The factorization scale dependence (denoted fi F ) may be exhibited explicitly by combining 
the second term in Eq. (|54]) with the corresponding term in ([43|): 



d(JAB- 



dQ 2 dQ 2 T 



V <x (H) 



ab 



N 



X(j)a/A(N + 1, fi F )4>b/B(N + 1,{1 F )t 

\E^(N,b,Q,Q) 
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r e ibQT exp 



=a,b 



exp 



D a f\N,b) 



(56) 



where is the full (lnA^ and constant) term in the iVth moment of the diagonal splitting 
function for parton i. Notice that we have set fi = Q in E^ eik \ In the resummed cross section (|5T)|), 

" 0^ 



dependence on the factorization scale \i F is suppressed compared to NLO [|5|, |26|, [31], \ 
because the eikonal anomalous dimensions compensate the evolution of parton distributions to 
leading power in N, at all orders in a s . 



Eq. (|56"D is our most general form in this paper for the Drell-Yan cross section, which we will 
approximate to NLL in Sec. 5. Applications to other processes are, as we shall see, conveniently 
carried out in terms of the coefficient function Q/j, Eq. (|34l). By comparing Eqs. (H) and (g3|), and 
once again using the explicit fi F dependence of the resummed exponent, we derive an analogous 
expression for the product of coefficient functions c^: 



■z a /a(N, b,Q,n) c b/b (N, h,Q,/j) = exp 



Q 2 du' 2 
//4 V' 2 % = a ,b 

i(eik)/ AT u n n\ i n( eik ) 



xexp E a r(N,b,Q,Q) + D^'(N,b) 



(57) 



This result will be useful for direct photon and other hard-scattering processes with factoring 
initial-state interactions, but with final-state color flow. 
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4. Single-particle Inclusive Cross Section 



In the following, we apply the methods outlined above to p^daAB-*j+x( x T) /dpr, the prompt 
photon inclusive cross section at measured pt- Unlike electroweak annihilation, however, we will 
not take the limit pt — * 0. As a result, we have no explicit logarithms of pt to resum in the hard- 
scattering cross section. Threshold resummation has been carried out for this process in [EE 



and its consequences studied phenomenologically in |26] , fig] , but there has been no generally 
accepted method to incorporate the kind of recoil corrections that are so important in describing 
the Iow-Qt limit of electroweak annihilation. These effects must, however, be present in the 
hard-scattering functions for this process at some level, even if they cancel almost completely. 
Our goal here is to develop a framework in which we can identify them systematically. This is a 
prerequisite to any reliable estimates of their influence. 

Models of "intrinsic" transverse momentum seem to suggest that recoil effects of magnitude 
similar to those in electroweak annihilation may be important |12, Recent experimental 



results appear difficult to explain without them JTTJ. The inputs to these analyses are primarily 



nonperturbative, in the form of Gaussian distributions in partonic transverse momenta, whose 
widths may be compared directly to nonperturbative parameters in electroweak annihilation [17 . 
At the same time, the interpretation of the theory and experiment is not without controversy 
15 , 16fl . For these reasons as well, it is of interest to reanalyze the prompt photon cross section in 



the light of the joint resummation [|T8j procedure introduced for electroweak annihilation above, 
in the formalism of collinear factorization. 



4.1. Hard-scattering functions for inclusive prompt photons 

The prompt photon cross section may be written in factorized form as 



p^da AB ^ y+x (x T ) r r da ab ^ 7 (x^,pr/fi, at a (ji)) 

j = dx a (Pa/A{X a ,fl) / dx b <p h / B (X b , fl) p T 

apT ab J J dpr 



where we define xt and x t by 



2 _ 4 PT ~2 _ A Pt /kqN 



S s 

with s = x a XbS. To leading power in N, moments oip\dajdpT are 

3 da ah ^ (N,p T /fx,a s (fx)) f 1 2 a2 n _ x 3 d<7 o6 ^ 7 (x^,p T /[M, a e (ji)) 



dpr Jo ' ' ' dpr 

= -~ L C dx\ ( x 2 T ) N -y T ^i. (eo) 

<j) a/a (N + 1, fi) c/> b/b (N + 1, /i) Jo TJ FT dp T 

Following our discussion of electroweak annihilation, we will use a combination of threshold and 
Qt resummation to estimate higher-order corrections in the partonic hard-scattering function, 
PT&ab-t-y/dpT, always working to leading power in the moment variable N. 

At fixed photon rapidity fj in the partonic center-of-mass system, singularities arise when the 
partonic center-of-mass energy y/B, reaches the threshold value of the final state necessary to 
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produce the photon at that pr and 17. Any such final state is kinematically equivalent to the 
photon recoiling against a massless jet. The minimum invariant mass at measured pt is 

4p 2 

Smin = ^Pt COsh 2 fj = . (61) 

This kinematic relation changes, however, when we take into account the recoil of the photon- 
jet pair transverse momentum, denoted Qt below, against perturbative initial-state radiation. 
In effect, when the incoming partons get a kick in the direction of the observed photon through 
initial-state radiation, the minimum invariant mass necessary to produce the photon at measured 
Pt decreases. These effects are present in higher orders in a, but are not reflected directly in logs 
of 1 — x 2 T . 

We can only begin to take recoil into account systematically, however, when we have defined 
what we mean by initial-state radiation, and what by the short-distance scattering. This is pre- 
cisely an issue of refactorization, discussed above in Sec. 2. There we showed that refactorization 
into generalized initial-state parton distributions, soft radiation and a short- distance process is 
quite natural to leading power in threshold singularities, (1 — z)' 1 , for electroweak annihilation. 
As we shall see, leading power at threshold for prompt photons allows refactorization in an anal- 
ogous manner, including now a new function for the final-state jet, and a generalization of the 
function for soft radiation. The essential point, however, is that at partonic threshold it is pos- 
sible to identify a 2 — > 2 short- distance scattering, involving only lines off-shell by 0(p T ), that 
underlies the production of the high-p<r photon. Our aim is to treat this short-distance scattering 
in the same manner as the short- distance cross section in electroweak annihilation, in terms of 
its invariant mass squared, s = Q 2 , and its transverse momentum (defined in the hadronic 
center-of-mass system). When is in the direction of the observed photon, s may be less than 
s, the partonic invariant mass squared in Eq. Q5S|). 

To take recoil into account in joint resummation, we identify the transverse momentum of 
the photon in the center-of-mass system of the final-state photon-jet pair as one-half the relative 
transverse momentum of the photon and recoiling jet at partonic threshold, 

Pt = Pt - ^ Qt • (62) 

In terms of the kinematics of the short-distance 2^2 subprocess, we can define a natural 
"scaling" variable for this refactorized scattering, analogous to xt and xt in Eq. (|59|) , 

~2 _ MPT? 1 (M , 

Q z cosh T] 

with fj the rapidity in the refactorized hard-scattering center-of-mass system. The variables x T 
and x T are related by 

4 = 4 (I ffi ' (64) 

a relation that we will use below in the analysis of moments. 

We will estimate the partonic cross section as the integral over Qt and Q 2 = s of the doubly- 
resummed cross section at measured Qt and Q 2 , limiting the Qt integral to a cut-off scale p. 
That is, we will write the cross section as 

, i (resum) 

"^^^ + A( PT , P) , (65) 



22 



with 

, (resum) , (resum) 

*^=J****«r*%*;W-™- (66) 

We have introduced a cut-off ft in the recoil transverse momentum Q^, which we include to avoid 
going outside the range where the approximations for joint resummation fail, that is, where the 
recoil transverse momentum becomes competitive with the hard scattering. The need for a 
matching condition for the resummed to fixed-order expressions at high recoil is familiar from 
QT-resummation in electroweak annihilation. Nonsingular finite-order terms, corrected for the 
matching P, 20 are included in A(pT,p)- The implementation of such a matching procedure 



remains to be carried out in the new joint resummation formalism. We shall, however, exhibit 
the theta function in in each of our expressions below, as a reminder of its importance. 

In our analysis below we will determine the jointly resummed cross section at fixed Q 2 and 
Qt- To develop a jointly resummed cross section, we begin, as for electroweak annihilation, with 
a study of perturbation theory near partonic threshold. 



4.2. Leading regions in single-photon production 



The derivation of a refactorization formula for the single-photon inclusive cross section is quite 
similar to the analysis that leads to Eq. (§) for electroweak annihilation. In the photon cross 
section, however, final-state interactions play an important role. To analyze their contributions, 
we need an analysis of the leading regions in the momentum space of cut diagrams that produce 
logarithmic corrections to the cross section. The analysis is quite similar to that carried out 
for threshold resummation in heavy quark |30J, |31| and jet production |3H}| , but now taking into 
account transverse momenta near threshold. In this analysis, we shall neglect, for the time being, 
fragmentation components in the prompt photon cross section, which are relatively modest in a 
significant kinematic region [p6 



The relevant leading regions for the 2 — ^ 2 partonic subprocess a + b — > c + 7 are illustrated 
by Fig. |2], which may be compared to Fig. [I] for electroweak annihilation. In addition to the 
jets 7Z a /a and TZb/b associated with the incoming partons, and the short- distance subdiagrams H, 
there is also a subdiagram J c that account for partons collinear to the outgoing parton c, and 
a new soft subdiagram S a bc, which accounts for soft radiation from the final as well as initial 
hard partons. These leading regions are of the general class discussed in Ref . [|J , and identified 
for hard-scattering cross sections in Ref. 
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on the basis of analyticity properties and power 
counting bounds. In the general case, there may be an arbitrary number of noncollinear jets in 
the final state. Leading power in the threshold variable, 1 — x 2 -,, limits their number to a single 
jet recoiling against the hard photon. 

The role of the final state jet J c in threshold resummation is well-understood, and has been 
treated in Refs. |32|, |33|, |26f . The outgoing jet may interact with soft radiation, through subdia- 



gram S a b c in Fig. 0. We may think of the soft function as associated with coherent soft radiation, 
describing emission and absorption by sources with specified four-velocities and color charges. 
The role of soft gluon functions encountered in threshold resummation has been extensively 
studied elsewhere [plj. We will give a formal definition of the fully factorized soft function in 



subsection 0~3 



There is a potential complication for the joint resummation of the single-photon cross section 
in the influence of soft radiation on the recoil transverse momentum of the photon- J c pair. By 
analogy to electroweak annihilation, we seek to resum logarithms of 1 — x\ as well as logarithms of 
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a 



b 




Figure 2: Leading region for prompt photon production near partonic threshold. The vertical 
line represents the final state, including the photon, 7, and the recoil jet, J c . 

the total transverse momentum of the partons involved in the underlying 2^2 hard scattering, 
Qt- This transverse momentum is defined relative to the (initial-state) hadronic center-of-mass 
system, evaluated at fixed observed photon transverse momentum, p?, in that frame. Referring 
to Fig. 0, we must ask, however, whether Qt is the same on both sides of the cut, i.e., for the 
amplitude and for its complex conjugate. An imbalance in the relative transverse momenta of 
the two hard-scatterings would make it necessary to introduce a more complicated convolution 
than in Eq. @. Such an imbalance, however, can arise only from the transverse momentum that 
flows through the soft function S a f, c between the incoming jet sub diagrams Ttf/f and the outgoing 
recoil jet J c . In the absence of such a flow, the transverse momentum of the final-state jet is a 
dependent quantity, and does not appear directly in the transverse momentum logarithms. More 
importantly, a single Qx describes the relative transverse momenta of the hard scattering on 
both sides of the cut. 

We must therefore study the flow of transverse momentum through the soft subdiagram. 
We begin with diagrams that have initial-state interactions only, in which soft gluons, as in elec- 
troweak annihilation, couple only to the initial-state jets, 1Z a / a and Hb/b in Fig. |2|. The transverse 
momentum that they carry into the final state must come from the hard scattering, through the 
initial-state partons, a and b, equally in both the amplitude and the complex conjugate ampli- 
tude, that is, with the same Qt on the two sides of the cut diagram of Fig. ||. At the same 
time, infrared divergences associated with coherent soft gluons cancel as we sum over different 
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(a) 




Figure 3: Diagrams that illustrate the cancellation of infrared divergences: (a) initial-state inter- 
actions, which require different diagrams, with different parton momenta k a and k^. (b) initial- 
final-state interference, for which the cancellation proceeds through cuts of a single diagram, at 
fixed k a and k b . 



cut diagrams, with different connections of soft gluons to the jet diagrams. In the diagrams 
necessary to cancel infrared divergences, therefore, the initial-state parton transverse momenta 
k a and k& will each vary. Examples are shown in Fig. |3]a. This imperfect match means that the 
cancellation proceeds through plus distributions in k^ and k^, and can produce logarithms of b 
in the impact parameter space conjugate to Qt- These purely initial-state interactions can be 
treated exactly as in Sec. 2 above. 

Consider next soft connections to the final-state jet, including interference between initial- 
and final-state interactions, as in Fig. [|b, where the soft function S a b c now has contributions in 
which soft radiation is emitted by an initial state line, and absorbed by the final state jet. As 
usual, we must sum over final states to cancel the infrared divergences associated with this soft 
radiation. There is, however, a crucial difference between the cancellation of initial-state and 
initial-final coherent soft radiation. In the case of initial-state radiation, as just described, we 
must sum over different diagrams to effect the cancellation. Final-state, and initial-final-state 
interference divergences, however, cancel in the sum over cuts of individual uncut diagrams, 
as discussed in context of factorization proofs in ||50|1 . Li and Lai also observed that the /in- 
dependence due to final state interactions cancels in the context of prompt photon production 
|17| . Hence, unlike singularities that arise from purely initial-state interactions, the cancellation 
of final-state singularities can be effected at fixed transverse momentum for all lines in a cut 
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diagram, in particular, for fixed pair transverse momentum Qt of the initial-state partons a and 
b, and hence for the short-distance process in both the amplitude and its complex conjugate. In 
contrast to initial-state interactions, the transverse-momentum dependence of the final-state soft 
radiation cancels algebraically, rather than through plus distributions. 

It is worth noting that the above result requires that we resum logarithms of the transverse 
momentum at the short- distance scattering, rather than of the observed photon-jet pair in the 
final state. In the latter case, the transverse momentum at the short- distance scattering depends 
on which of the soft gluons attached to the final-state jet are virtual and which are real, and 
the cancellation of final state infrared divergences reverts from algebra [51J to plus distributions, 
and may produce logarithms in impact parameter space. When the recoil jet is observed inde- 
pendently, therefore, a somewhat different analysis is necessary, which we shall not carry out 
here. 

To summarize our considerations so far, only the transverse momenta of coherent soft radia- 
tion associated with initial-state hard partons must be taken into account in joint resummation. 
In contrast, final-state and coherent soft radiation linking the outgoing jet with one or both of 
the incoming jets produces no logarithms in the pair transverse momentum. 

The cancellation between the final states in Fig. still requires an integral over energy || . 
For this reason, both initial- and final-state interactions contribute logarithms to threshold resum- 
mation. We must incorporate the distinction between initial-state and initial-final interference 
logarithms into the refactorized convolution that generalizes Eq. @ for electroweak annihilation 
to the case of prompt photons. In the next subsection, we show that it is possible to do this by 
separating purely initial-state soft radiation from initial-final interference, at least to the level of 
next-to-leading logarithms in both threshold and pair transverse momenta. 



4.3. The soft function 

To derive an analog of the electroweak annihilation refactorization formula, Eq. (|Sp, for direct 
photon production, it is necessary to identify a function that summarizes final-state soft gluon 
radiation. In particular, we want to separate those effects associated with initial-state radiation, 
which are sensitive to transverse momenta, from those from the final state, which are not. In 
this subsection, we will construct such a function. 

The soft function will be built from nonabelian phase operators, following the discussion of 
Sec. 2.3 above. We begin by generalizing Eq. (|20|) , which describes the annihilation of a pair of 
phase operators to a product that corresponds to the color flow |J{], |3T], |ST|| at the hard scattering, 
for the partonic process, a + b — > 7 + c: 



$g(oo,0;X) 



sg>(0,-oo;X) 



°abc )d c ;d b ,d a 



d b ,e b 



da )6 a 



(67) 



As in Sec. 2.3, we go on construct an eikonal cross section in the form of Eq. (pl|), at fixed 
soft-gluon energy, parameterized as w s Q, and fixed transverse momentum k [[3(1 EM, 
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(68) 
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where the trace is over the external color indices (e*) of the operators. In this expression, h is 
the unit vector in the time direction, and we leave the renormalization scale \x free. The relevant 
transforms with respect to w s and k are Laplace and Fourier, respectively, 

*2? f^T> bQ, e)= r dw s e- Nu ' ! rf 2 ke- b ' k a^\w s , Q, k, e) . (69) 



Q 

The soft radiation function for the a + b — > 7 + c subprocess in prompt photon production is now 
constructed from <j( eik ) ( N, b) , by analogy to the soft function U for electroweak annihilation, Eq. 

We find the soft function by dividing the transformed eikonal cross section (|69|) by functions 
that eliminate double counting with the external jets near partonic threshold: both incoming, 
and, in this case, outgoing. For the incoming lines, these functions are the TZ^ 1 ^ ( N, b) , defined in 
Eq. (|30|). Similarly, for the outgoing jet, we identify a new partonic function, which will appear in 
the refactorization formula, along with its eikonal partner. As we have seen, infrared divergences 
associated with final state interactions cancel at fixed recoil for the hard scattering. We may 
therefore treat the outgoing jet inclusively in its transverse momentum. The relevant functions 



are then the same as those encountered in pure threshold resummation ||34j| . For example, for a 
quark jet they may be defined as two-point functions, 

J^ cik \Nfi/Q) = fdw c e~ Nwc f da e- i{ ^ Q/V ^ )cT 



x — Tr(0|^t(0) s£V&)|o> 



J c (Nfi/Q,Q/fi) = Jdw c e- Nw < j dXda 



" Q Tr(0|g c (0) 7 • & g c (A& + a(3 c + y)|0) 



V2N C 

dw c e- Nw < J c (w c ,Q) , (70) 

for the eikonal and partonic jet, respectively. In the partonic jet function, j3 c is the light-like 
velocity vector in the direction of the jet, and (3 C is the velocity vector opposite to (3 C in the 
overall partonic center-of-mass. The matrix elements are evaluated in n ■ A = gauge. The 
traces refer to color and Dirac indices. Outgoing gluon jets may be defined analogously, with 
operators similar to those used for fragmentation functions | ^9[| . In momentum space, w c Q 2 is 
the squared invariant mass of the outgoing jet. 

The above reasoning leads to the following generalization of the soft function to prompt 
photon processes: 



g(ab^c) l_Q_ h Q a (u) n ) = (^,b,a a (/x), v 

J n^ ik \N l i/Q,bQ,e)nt ;ik \N l i/Q,bQ,e).n cik \N t i/Q)' 



As in the case of electroweak annihilation, Eq. (p7|) , this expression is a rewriting of the refac- 
torization of the eikonal cross section into incoming and outgoing jets and soft radiation. Again, 
the K K : r remove initial-state eikonal radiation from <j\J^ , and J c ( eik )(iV) removes dependence 
on the outgoing jet. In this form, however, the soft function inherits the gauge dependence im- 
plicit in Eq. ([H]) of the incoming and outgoing jet functions. We can eliminate this dependence, 



27 



and simplify the overall formalism, by following an observation made in Refs. |3(| [51]] and em- 
ployed in |j32| . We define a variant soft function in transform space by dividing S by a factor 

— 1 /2 — 1 /2 

Ug (Nfji/Q, bQ) for each of the incoming jets a and b, and by Uj (N, 0) for the outgoing jet, 



g(a^c) (_fi_ >b Q,a fl (//),7i 



UlL\Nn/Q, bQ) U^iNfi/Q, bQ) uU^N^/Q, 0) ' 



(72) 



By shifting factors of U 1 ^ 2 from the soft radiation function to the jets in transform space, we 
produce slightly modified jet functions, of the convolution form 

1Z' a / a (x, Q,k,e) = J dwdy 8(w + y — x) J dk r dk u 5 2 (k r + k u + k) 

1 /2 

x Ti a / a (w,Q,k r ,e)U a i (y,Q,k u ) 
= [ di o/a (£, ») ~ c a/a (|, k, Q, f?j + 0(1/N) , 
j' c (w c ,Q) = J dw'dw u 5(w' + w u -w c ) J c (w',Q) J dk u U^ 2 (w u , k u ) , (73) 

where in the second equality for lZ' a / a , c a / a is the double inverse transform of the infrared safe 
function c a / a defined in Eq. (|34|) above. 



4.4. Refactorization, recoil and the resummed cross section 

The refactorization formula for prompt photon production generalizes the corresponding expres- 
sion for electroweak annihilation, Eq. (|9|), by including the outgoing jet function J' c in Eq. (|73|), 
and the modified soft radiation function 5", Eq. (0). In transform space the refactorization is 
in terms of products; in momentum space in terms of convolutions, 

dcr (re ^ um) f f 

P ^ dQ 2 (PQ 1 dp = J dx * d2k * n 'a/a(xa,k a ,Q) J dx b d 2 k b K' b/b (x b ,k b ,Q) 

x J dw c J' c (w c , Q) J dw s J d 2 k s S'H-tc) Q ? ks; as(/i)) 

x i 6(1 - Q 2 /S - (1 - x a ) - (1 - x b ) -w s - w c ) 5 2 (Q T + k a + k 6 ) 

., riab^c, , \ ~2 \ 3 dp' T da (Xt) 

XC S -(a s (^x T )p T — ^ . (74) 

Eq. ([74] ) generates, order-by-order in perturbation theory, the same singularities at Qt = as 
the full cross section, to leading power in 1 — Q 2 /S. The factor dp' T /dpT compensates for the 
difference in phase space between fixing pt and pt ■ The function Q^ h ~*^ c ^ — \ _|_ 0(a s ), times 
the Born cross section, is the perturbative short- distance function, which in the case of prompt 
photon production (as opposed to electroweak annihilation) depends on x\. (] This means that, 

5 From Eq. (||), Xrp IS determined by pr, Qt and Q. 
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beyond the lowest order, the short- distance scattering function need not have the same angular 
dependence as the Born cross section. The short-distance function contains only corrections 
from (virtual) lines that are off-shell by at least 0(pt)- It contains no real-gluon emission, since 
all radiation has, to leading power in 1 — Q 2 / S, been absorbed into the long-distance functions 
1Z', J' c and 5". The computation of these short-distance functions is equivalent to the matching 
conditions of effective field theories. 

With S' constructed as above, the soft transverse momentum k s is associated entirely with 
final-state interactions, and is not included in the recoil momentum of the hard subprocess. We 
may therefore integrate over k s and redefine 



where on the left we retain the same notation for the function, but omit the transverse argument. 
Our refactorization formula then simplifies to 

Pt d Q 2 ^2Q^ pT = J dx a d 2 k a Tl' a/a (x a , k a , Q) J dx b d 2 k b TZ' b/b {x b , k b , Q) 
x / dw c J' c (w c ,Q) [dw s S'( ab ^(w s ,Q,a s (ri) 



x i 5(1 - Q 2 /S - (1 - x a ) - (1 - x b ) -w s - w c ) 
x 5 2 (Qr + k a + k 6 ) 

w n {ab^ c ), / > ~2\ 3 dp' T da^ (x T ) 

xQ Mfi),x T )p T - -ri . 76 

dpx dp T 

To simplify the notation further, we also introduce a new function J 7 , which combines contribu- 
tions from the final-state jet and soft final-state radiation, 

Fabciwf) = / dw s dw c s"( a6 -»Tc) ( Ws) g ; a s (/i)) J' c (w c , Q) 5(w f - w s - w c ) . (77) 



As described above, this convolution does not involve the recoil transverse momentum. In mo- 
ment space, Eq. ( fT?D becomes a product 

FaUN) = S' (o6 ^ c) (Jfy J' C {N) + 0(1/N) . (78) 

We will discuss the explicit ^-dependence from soft gluon radiation and the final-state jet below. 

Eq. (|76|) defines recoil in the prompt photon cross section in much the same way that Eq. 
(|]) defines it for electroweak annihilation. The short distance functions, C$ dcr^ ' (x T ) / dpx in- 
clude only lines off-shell by 0(p' T 2 ), and are evaluated at zero relative transverse momenta for 
initial-state partons. Expansions in the transverse momenta of the incoming partons are to be 
absorbed into higher orders in Cs- In this convolution form, however, the kinematics of the hard 
scattering influences the cancellation of singularities at vanishing k a and k^,. This procedure 
has a straightforward interpretation order-by-order. At fixed order, all contributions that are 
singular at threshold for fixed Q 2 may be put in the form of Eq. (76). To evaluate the cross 
section at any fixed order, we would integrate each such contribution over Q 2 and Qt, as in 
Eq. (|66D , with no further approximations. The result would contain finite corrections resulting 
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from the kinematic linkage of the hard scattering with the cancellation of singularities in trans- 
verse momentum. We would then sum to all orders. In the resummed formalism, we simply 
approximate the short-distance function to 0(a s ), sum the singularities to all orders at leading 
or next-to-leading logarithmic accuracy, and do the Q 2 and Qt integrals last. In the following, 
we shall derive the consequences of this reorganization of perturbation theory. 

We are now ready to return to Eq. (|66D , and derive the jointly-resummed partonic cross 
section as an integral of the differential resummed cross section ( |76|) over the hard-scattering 
scale Q and its relative transverse momentum, Q^. Changing variables from Q 2 to we derive 

■3 ; (resum) a 

Prdcr K ab _^ p^ 



dpT 8nS 2 J2 J Q d ^ d ^3 ti/aiZh^j/biZj^) 

J Jo . h _ ^2 



J <i 2 kj d 2 ~kj S 2 (Qt + kj + kj) J dxi J dxj J dxf 5 



~2 [ Pt 2 \ X T 

T ^2 



p\ J XiXj(l — Xf) 



( ' Xi 2p' T \ _ ( Xj , 2p' T \ . ~2\ { 1 \ 

xdi/i W "^T J/3 W ^ x " Xt) { x t x 3 {i-x f) ) ' (79) 



where we have used Eq. ( |73D to isolate the partonic hard- scattering function in MS scheme, 
and Eq. ([77D to summarize the contribution of soft and jet radiation in the final state. In the 
argument of the delta function, we have used Eq. ( |64D to reexpress the ratio Q 2 /S as 



S \ Pjr 1 J 



(8°) 



In Eq. (|7Up, we have also reexpressed the Born cross sections in terms of the 2^2 invariant 
amplitudes M(x T ), and have used the approximate relation p' T (dp' T / dpx) = Pt, valid up to a 
nonsingular term of 0(Qt), which we neglect. 

It is a relatively small step from Eq. ( [7§D to a jointly resummed cross section for hadrons, by 
replacing partonic by hadronic MS distributions, 

„3 J. (resum) 4 

p T aa AB ^ _ p T 



dpT 8nS 2 12 J d & d £j < Pi/A(£i,V)<Pj/B(£j,V) 

K / 7T^¥ Qt / dx T f C\ \a s (fi),x 



f Xi 2p' T \ _ f Xj 2p' r \ _ 2 / 1 \ 

x di/i w ' ~ Cj/j w ' ^1 -xj) ) ' (81) 

where we have also replaced the convolution in transverse momenta by a Fourier transform, so 
that the functions c are now in impact parameter space. 

Eq. (|8lD factorizes under x\ moments at large N, up to 1/N corrections. Thus, following 
essentially the same steps that led to Eq. (|43|) for electroweak annihilation, we find the physical 
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cross section as an inverse transform Hl9 



pWabT] P 4 t ^ fdN~ 



dpT 8nS 2 Jc 2m 



- ° v 7 \JI-Xt 
f d 2 Q T , f S \ N+1 

x J d 2 b e lh QT Ci/i (n, b, ^ c j/:i (n, b 



// F ijk {N,x' T ). (82) 



In the c's, the impact parameter b is conjugate to Q T at fixed values of Q 2 = 4p' T 2 /x^. In 
this form, however, each Cj/j also depends implicitly on Q^, through p' T 2 . Explicit forms for 
the functions Q/j have been given in Eq. ( |57D above. In the next subsection, we will find the 
iV-dependence of the final-state function T. 



4.5. Resummation for the final state 



Constructed as in Eq. (0), £'(<^7c) sa ti s fies the same renormalization group equation as the 
soft functions encountered in heavy quark [RH] and jet [R4| threshold resummations, 



dji 



a s (jj) 



-2ReT 



(afe-»7c) 
S' 1 



13) 



where the anomalous dimension r^ 6 ^^) j s a function of the velocities /3j associated with the phase 
lines corresponding to partons a, b and c. For the special case of prompt photon production, Tg' 
is a number, rather than a matrix, because the short-distance cross section has only one color 
The solution to (|83|) is therefore a simple exponential, 
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C</(a.b— >7c) 



S'^^(l,a s (Q/N)) 



x exp 



—2 Re r v 9 , 
/' /'' 



M) 



The logarithm of the soft function, constructed in this fashion, has at most a single logarithm of 
N per loop, so that by calculating p^ b ~ K >' c ) to one loop, we determine S'^ ab ~*' yc ^ at the level of its 
leading logarithms, that is, ln L N in the exponent, while remaining at NLL in the overall cross 



(ab— »7c) 
S' 



is calculated from the MS counterterms for the three diagrams shown in Fig. 



section. T 

in which a single gluon is exchanged between pairs of eikonals. The details of these calculations, 
which are to be carried out in an axial n ■ A = gauge, are described in |5T] . Many terms cancel, 
and in Fig. [|a, the additional crosses on the gluon lines denote a slightly modified A = gauge 
propagator: 

1 / n^k u + k^n u 1 k>*k 



k 2 



n ■ k 



in 



kf 



5) 



31 





(a) 




(b) 



Figure 4: (a) Diagrams which contribute to the soft anomalous dimension. The crosses on the 
gluon lines denote the modified propagator given in the text. The remainder of the A = gauge 
propagator cancels the real part of (b). 



The factor 1/2 in this expression is specific to the lowest-order calculation. It takes into account 
the effect of subtracting the factor U aSi U b i implicit in the definition of c a q,, Eq. ([34]). The "missing" 
terms in Eq. fl85|) completely cancel the real part of the diagram of Fig. |]b. The resulting 
anomalous dimension is then 



ReT 



ReT 



ill- 
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2^ 



2C F In (-j- 
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16) 



where s, t and u are the invariants for the partonic 2 — > 2 subprocess. These anomalous dimen- 
sions are independent of the gauge vector n M at one loop [[51]]. The functions S'^ ab ^" ,c \ defined 
as the solutions to Eq. fl83|), with the usual boundary condition S'^^^^N = 0) = 1, are free of 
initial-state radiation that would be sensitive to recoil at the logarithmic level. They continue to 
contribute to the threshold phase space through the energy w s Q []3(| j34"l , in a manner described 
below. 

Explicit expressions for S'(N) are found from Eqs. ([84] ) and (|86| ) above. The transform of 
the recoil jet function may be found from Ref. [52]|, and is given in its most general form by 



J' C (N) 
E' C (N) 



exp[E' c (N)' 



JV-l 



dz- 



dX 



1 - z 



A c [a s (V\Q)] 
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B' c [a s (VT^Q)] + B" c [a s ((l - z)Q)\) , 



(87) 



where the A c are given in Eq. ([T5|). In the specific normalization chosen for the recoil jet, J' c , in 
Eq. ( f73l) above, .B" vanishes, while 




These results, along with Eq. (|57|) for the functions Cj/j, specify the explicit NLL resummed cross 
section as an inverse Mellin moment in Eq. (|S2|). We will give the relevant expressions in Sec. 
5. We close this section with a discussion of the additional considerations necessary to include 
fragmentation effects in the formalism. 



4.6. Including fragmentation 

For single-particle inclusive hadron production, and even for high-p^ photon production, we must 
supplement the considerations above to include final-state fragmentation. This is relatively easy, 
since the dynamics of fragmentation factorizes from initial-state, hard virtual, and soft coherent 
radiation as well as from the final-state corrections associated with the recoiling jet, as illustrated 
in Fig. |5|. The analysis of this section applies to pure threshold resummation (as discussed in 
Pfifl) as well as to joint resummation. 

The underlying short- distance scattering subprocess is again a 2 — > 2 reaction, but now with 
two outgoing colored particles, one of which fragments into the observed hadron, which we may 
take to be a pion (or photon): 

a(x aPa ) + b(x bPb ) -> c + d(P/z) -> tt(P) + X , (89) 

where we have exhibited the partons' momenta. The outgoing jet, Jd, fragments into the observed 
hadron (jr). The short-distance scattering involves more than one color flow, the set depending 
on the flavors of the partons in Eq. (j55|). At any order in perturbation theory, the color flow in 
the amplitude and complex conjugate need not be the same. Referring to Fig. [5], and adopting 
the notation of electroweak annihilation, we represent the (dimensionless) short- distance function 
as pip dajj(x^)/dp' T = h*j hi, where f denotes the 2^2 partonic scattering of Eq. (|89f) . The soft 
function is built in the same way as for prompt photon production, treating the color tensors ci t j 
for each color exchange at the hard scattering as an effective local vertex linking phase operators 
in the flavors of the partons a. . .d. The soft function is written as Sjj . The same arguments 
regarding the cancellation of recoil effects in the soft function apply to single-hadron inclusive 
as to prompt photon production. Also, the normalization of the soft function follows Eq. (72") 



above, for each of the matrix elements of SjJ, but now with factors U 1 ^ 2 for both outgoing jets. 
The resulting anomalous dimensions for most relevant flavor flows have been computed in Refs. 
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Near threshold, in the notation of Eqs. ( [741) and ( j30|) , the phase space is slightly modified, 
relative to the prompt photon case, because the fragmenting jet has invariant mass gj > 0, and 
because the outgoing parton at the hard scattering carries momentum P/z > P, which shifts the 
scaling variable xt by a factor z. This changes the transverse momentum in the hard-scattering 
center-of-mass system: p' T — ► p'y(-z) = Pt/z — Qr/2, which is related to the short-distance 
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Figure 5: Leading region for a single-particle inclusive cross section in cut diagram form, cj and 
Cj represent color tensors |)T . 

scattering and overall invariant masses squared, Q 2 and S, by 



In the second relation we have expanded in 1 — z, and have neglected corrections of the form 
(1 — z)Qt- 111 place of the argument of the first delta function in Eq. fl74}) for resummed prompt 
photon production, we find the following kinematic constraint near partonic threshold for single- 
particle inclusive cross sections with fragmentation: 

O 2 

1 = y + (1 - x a ) + (1 - x b ) + w s + q 2 /S + q 2 JS 

1 2 2 

" " + 2(1 - z) + (1 - x a ) + (1 - x b ) +w s + q 2 /S + q 2 /S , (91) 





2 ™2 


\Pt 





with x| as in Eq. ( |59|) and £^ as in (|63|) . 

The refactorization that generalizes Eq. ( |7^ ) includes an additional integral over z, linked 
through the restriction (|9l|) , and a function that describes the fragmentation of parton d into 
hadron 7r(P). We parameterize the momentum of parton d as 

S^^^+p^S ( 92 ) 

where ^ is the opposite-directed unit light-cone vector relative to momentum P M = (P • fid) (3%, 
which is itself taken to be lightlike. The second form introduces a dimensionless variable Wd = 
ql/S. 

The fragmentation dynamics of the outgoing jet can be described by a function that is quite 




similar to a standard fragmentation function [29|, and to the inclusive jet functions (fZOQ above. 
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As above, we illustrate the case of a quark qf, computed in n ■ A = gauge, 



Q r dXda 



f/n ' ' ' ' 8N C Y J 271 



c 



-i(P-0d/z)>^-iqd-0dO- 



x Tr <0| q f (0) \f(P),X) fa (f(P),X\q f (\/3 d + afa) |0> 

f 1 d£ ( z \ 
= J 76/// 7,%^? (93) 

with //pa factorization scale. In the first, defining, relation, the trace is over both color and 
Dirac indices. A sum over spins is understood. In the second equality, -D/// is a fragmentation 
function, which we may define in MS factorization scheme. As usual, at partonic threshold, 
the infrared-safe coefficient function ed/d may be taken diagonal in flavor, up to corrections that 
vanish as 1/N in moment space. 

It is convenient to define a final state threshold function by analogy to Eq. flTT|), as a convo- 
lution of the soft function with the recoiling jet, 



J 1 



{Wf) ^\l dWs ^ 6 ~ WS ~ q2jS ) S ' Jlii) K ' Q ' J ^e) • ( 94 ) 

Laplace moments reduce the convolution to a product, 



J 



= ^ (f) (]^'^^)) J'^)- ( 95 ) 

It is now straightforward to generalize each of Eqs. ([76]), fl82|) and (|8lD , corresponding to refac- 
torization at partonic threshold, and to the jointly-resummed hadronic cross section written as 
an inverse moment, and in convolution form, respectively. We give below the first two of these 
generalized expressions explicitly. 

The refactorized partonic single-particle inclusive cross section is given by 

dai Tesnm) 



P ^ dQ 2 °d 2 Q T dpr = I dXadj2 ^ a n 'a/a( x a,K,Q) J dx b d 2 k b n' b/b (x b ,k b ,Q) 

x 5 2 (Q T + k a + k fe ) J dz J, 

E /*»/*#(»/> rf^ dp . 
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dw d 
dp' T da { j](x 2 T ) 



x 



r 1 d£ ( z \ 
J -7- e d / d I -,w d ,P,ii F \ D d / d {£,n Fl e) 



x i 5 (l - Q 2 /S - (1 - x a ) - (1 - x b ) - Wf - w d ) , (96) 

where xt is specified by Eq. (p0|). As indicated in Eq. (|66|) , the resummed cross section is found by 
integrating over Q 2 and Qt- The form convenient for a; T -moments is found by changing variables 
from Q 2 to x\. We compute the hard-scattering cross section in moment space, by dividing by 
perturbative distributions for the incoming partons, and by a perturbative fragmentation function 
for the parton that fragments into the observed particle. 
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The resummed hadronic cross section is found by replacing the perturbative distributions and 
fragmentation function by their hadronic analogs in moment space, and inverting the transform 
with respect to x\. The result is 

„3 J n (resum) 4 



x 



(2^ 9 "-'^> Up r Q r /2p J 2 (97) 
x | rf 2 b e ib ' Qr c a/a ( JV, b, M J c b/b hv, b, jF<?(iV, 4) e d/d (2iV + 3, N, P, fi F ) , 

where we define 

e d/d (M, N, P, fi F ) = Jdy dw d e -m-y)-N Wd ^ ^ ^ p> ^ (gg) 

In Eq. (p8t), M)j denotes the projection of the amplitude squared at lowest nonvanishing order 

(f) 

on each of the color flows in the amplitude and complex conjugate, and C s is the corresponding 
expansion in a s of the complete hard-scattering function. A similar form was found in Ref. |26| 
in the context of threshold resummation for the fragmentation component of high-p^ photon 
production. The function e d u in Eq. (p8|) is a double moment, with respect to the jet invariant 
mass squared, w d Q 2 and the scaling variable, y. Its double logarithmic behavior [26], however, 
is determined entirely by the latter, because in the momentum configurations that give rise to 
double logarithms, we have 1 — y ^> w d , corresponding to the collinear emission of soft gluons. 
Finally, we note once again the dependence of the resummed cross section on a cutoff scale /L At 
this scale, the resummed cross section must be matched to a finite-order, or partially resummed 
cross section. The investigation of the best implementation of this procedure remains for future 
work. 



5. Exponents at NLL 



In this section we apply joint resummation to electroweak annihilation and prompt photon pro- 
duction, and exhibit the relevant exponents to next-to-leading logarithm. For prompt photons, 
these expressions were the basis of the phenomenological estimates given in Ref. |19[ . We reserve 



for future work the corresponding results for single-particle annihilation including fragmenta- 
tion, since they will require slightly more elaborate calculations involving color mixing at the 
hard scattering |3(], [51| . 



5.1. Electroweak annihilation 



Starting from Eq. (|56|) , we can identify an explicit expression for the resummed electroweak 
annihilation cross section that is accurate to NLL in both N and b 
the exponent D a b in (p6|) 



19 



We recall first that 
contributes only at NNLL. For the NLL exponent E a b, we have used 
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for guidance the NLL approximation to threshold resummation introduced in Ref. 
expression for the perturbative exponent E^ b T in Eq. fl56|) is 



E™(N,b,Q,Q) 



(b/a+N/Q) 



% [A a (a s (^)) + A b (a s (^))] 2ln^f 



Our 



(99) 



where as above, we define N = Ne lE . Eq. (|9~9| ) approaches the normal forms of QT-resummation 
(in 6-space) as b — > oo at fixed N, and of threshold resummation for N — > oo at fixed 6. We have 
introduced an explicit dimensionless scale c\ into the 6-dependence of the lower limit in the first 
integral. In Eq. ([99]) , the factorization scale, set to \x F = Q in the upper limit of the integral, 
replaces the conventional upper limit c<iQ in the Qy-resummation formalism || [|. 

Following the format of | 26| , we find the following closed expression for E^, accurate to NLL 
in both iV and b. 



E£(N,b,Q,fM F ) 
where for a = q or g we define 

and 
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(101) 
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+ ln(l - 2/3) 



(102) 



in terms of moment variables, 



A = b a 8 (ji) In (Nj , 

/? = 6 a s (/i) In (N + bQ/ Cl ) . 



(103) 



The coefficients A® and y4^ 2) in Eqs. ( jlffi] ) and ([TU2] ) have been given in Eq. (^3j). The last 
term on the right-hand side of Eq. ( |102j ) is the contribution at NLL of the first exponent of Eq. 
(|56|), including the anomalous dimensions. The beta function coefficients in these expressions are 
given by 



b 
bi 



IIC A - AT R N F 
12vr 

17C\ - 10C A T R N F - 6C F T R N F 
24vr 2 



(104) 



We recall that Eq. (|56| ) was derived choosing the renormalization scale as Q, the mass of the 
produced boson. In the expressions immediately above, fi F is the factorization scale, while the 
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argument fi denotes the renormalization scale dependence that arises when the scale of a s is 
shifted from Q to \i. 

Using these results in Eq. (f)(J), we derive doubly-resummed expressions for the production 
of W, Z and Higgs at measured Qt, through electroweak annihilation, accurate to NLL in both 
transforms. Compared to existing Qr-resummation formalisms [p0| , pl|j , we anticipate modest 
changes associated with the additional threshold resummation, especially for the W and Z, when 
their mass is far below collider energies, as at the Tevatron. Nevertheless, we expect a decreased 
sensitivity to the choice of factorization scale |25|, |26|, £|7|, [48, f49|. At the same time 
be of interest to study the 6-space integral in 




it may 
" prescription 



introduced in Ref. 



in the "minimal-principal value 
We shall not pursue these phenomenological implications here, however. 



In addition to a perturbative exponent, we expect that nonperturbative contributions, of the 
sort familiar from Qr-resummation, will be phenomenologically important whenever Q F <C Q 
HH, ^TJ. As noted above, the form of power corrections associated with the running coupling 
may be inferred from Eqs. (p3|)-(p5|), with the result that only even powers in Q~ l and b are 
required p2l |4 



5.2. Prompt photons 

We now turn to the prompt photon cross section, and show how the NLL resummed cross section, 
already studied in a preliminary fashion in Ref. JB5|, is derived. We start with Eq. (|5!?D, which 
expresses the resummed cross section in terms of doubly-transformed initial-state coefficient 
functions c, given in Eq. (f)7]), and final-state functions JF, specified by Eqs. (|78|) , (|84l) and (|37l). 
All logarithms of b and N exponentiate, and we find the form 
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(105) 



The final factor in this expression, the inverse transform of the exponential, was described as a 
"profile" function in Ref. jB|, and was denoted Pij{N, Q T ,Q,fi). 
The exponential moment dependence at NLL is given explicitly by 



E^ lk (N,b,Q,ti F ) = H t (N J b,Q jf i F ) + H J (N } b,Q,^ F ) + F k (N,Q) + G i:jk (N). (106) 

Here the contributions from the initial state jets in q/j are computed from Eq. (p9|), and are 
contained in the functions 



H a (N,b,Q,fx F ) 



a s (n) 



h^(X,(3) + h^(X,(3,Q, f ,, f , F ) 



(107) 
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with the h$ given above in Eq. ( |102| ). Contributions in from the final state jets are the LL 
and NLL functions 

F a (N,Q) = ^-/i°)(A) + /i 1 )(A,Q,/i), (108) 

while those in Tijk from the soft functions, are NLL only, 

G abc {N) = g^ c (X). (109) 

The new functions /M are found from Eqs. fl3T|) and (|38"D, while the g^} c are computed following 
Eqs. (|84f) and (|^), with the results: 

/f(A) = 2^°) (A/2, A/2) -4°) (A, A), (110) 
/«(A,Q,/i) = 2/£> (A/2, A/2, Q tf i,Q)-h^ (A, A, Q,/i,Q) 
/l (1) ln2 

+ _2_L(l I1 (l-2A)-lB(l-A))--^ln(l-A), (111) 

7T&0 TT&o 

<?£(A) = ln(l-2A)lnf^ , (112) 
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C F \n[— \ + ]-C A \n (- 



u 



(113) 



where the coefficients A$ and B^> have been defined in Eqs. (|T5| ) and (0), respectively. These 
formulas were employed in Ref. jnj to test the sensitivity of the inclusive prompt photon cross 
section to joint resummationQ. 

Again, to the perturbative expressions, nonperturbative corrections must in general be added. 
For prompt photons, we do not yet have all-orders expressions of the full generality of Eqs. (|53|)- 
fl55| ) for electroweak annihilation, so we cannot yet conclude that all nonperturbative corrections 
are of even powers in b and N. We consider this the most likely scenario, however, and in Ref. 



19 1 , the significant role of a modest nonperturbative term proportional to b in the exponent 



was highlighted. Clearly, this issue bears further study. 



5.3. Recoil and enhancement 

The jointly- resummed electroweak annihilation cross section, Eq. fl5E|), differs from 6-space for- 
malisms for Qr-resummation g g(J j n the inclusion of threshold logarithms at leading power 
in the Mellin moment variable N. At NLL, these corrections are included through the modi- 
fied lower limit in the exponent, Eq. (|99"1), which stabilizes the exponent to an integration range 
Q/N<fi'<Qa.tb = 0. This limit gives exactly the threshold resummed exponent for the 
(^T-integrated cross section |5|, 

In view of the above, we see that, in the \J integral in the exponent, Eq. (|9~9"D , the effects of 
QT-resummation and of threshold resummation are additive. This feature may be understood 
from the distinct origins of threshold and transverse momentum logarithms. The logarithms of 
iV in threshold resummation, which enhance the cross section, come primarily from subtracting 
negative corrections to the perturbative parton distributions in the construction of the hard- 
scattering function in the MS (or other) scheme ||. These subtractions produce enhancements 

6 We note a slight difference in notation from Ref. [[[9) , where the functions Cs were not exhibited separately. 
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because virtual corrections suppress the denominator of Eq. (£|) at the edge of phase space, i.e., at 
large N. The logarithms of b, however, are associated with real-gluon emission, in the numerator 
of Eq. (|]), the partonic cross section. They are negative in the exponent, because they cancel the 
divergent lower limit of the virtual // integral in Eq. fl99f) . Notice that the integrand changes sign 
at fi' = Q/N, which is the dividing line between threshold and transverse momentum dependence. 
Very roughly speaking, it is possible to resum threshold and transverse momentum logarithms 
simultaneously because they come from "different" gluons: gluons from the parton distributions 
for the former and from the unsubtracted partonic cross section for the latter. In essence, when 
recoil is taken into account, the Sudakov suppression associated with large values of b in ( P§D 
does not cancel the threshold enhancement, but redistributes it over a range of partonic invariant 
masses and momentum transfers at the hard scattering. 

In the jointly-resummed prompt photon cross section, Eq. (|105 ), the explicit enhancement of 
the integrand is associated with the inclusion of recoil, relative to the threshold-resummed cross 
section. The recoil allows x T < x\ for pt < Pt, even while the full partonic invariant s = £i£2*S' 
is bounded from below by 4p T . In moment form, Eq. fl82|) , this enhancement appears through 
the factor (S/Ap' T 2 ) N > (S/4p T ) N . Eq. ( 105 ) reduces to the inclusive threshold-resummed cross 
section in Ref. [26| if we neglect recoil, that is, if we set b to zero in the functions Cj/j, or at NLL, 
the exponent Eij in Eq. ( |106| ). Then, the exponent reverts to its form in threshold resummation, 
and the b integral gives a delta function, which eliminates the Qt integral, and freezes = pt- 



6. Conclusions 

We have described a joint resummation procedure for threshold and transverse momentum sin- 
gularities. The resummation organizes a well-defined set of corrections to single-particle and 
electroweak annihilation cross sections at measured transverse momentum, to all orders in per- 
turbation theory. Although the arguments for specific applications are somewhat involved, the 
basic observation is relatively simple: for those contributions that are singular at partonic thresh- 
old, the transverse momentum, Qt, of the short- distance scattering can be identified meaning- 
fully. The factorization properties of perturbation theory near partonic threshold allow us to 
control logarithms simultaneously in Qt and the relevant threshold variable, 1 — z, for elec- 
troweak annihilation, or 1 — x T for single-particle inclusive (1PI) cross sections. For electroweak 
annihilation cross sections, Qt may be identified with the observed final-state boson. For high-p^ 
1PI cross sections, Qt must be integrated to derive the hard-scattering function in the formal- 
ism of collinear factorization. The integral over singular distributions in Qt leaves behind finite 
remainders that modify, and may enhance, the predictions of threshold resummation. 

Joint resummation extends our control over a class of effects that can have an important 
phenomenological influence. At the same time, the resummed expressions afford new insights into 
nonperturbative power corrections in hard-scattering cross sections, which may be competitive 
with resummed perturbation theory, or even overshadow it. A preliminary study points to the 
importance of both high orders in perturbation theory and of nonperturbative effects in prompt 
photon production at fixed target energies [|H| . 

For electroweak annihilation at measured Qt <C Q, the gap between existing formalisms that 
resum in transverse momentum, and the joint resummation described here, should be relatively 
easy to close. Such an application would produce, we conjecture, a decrease in sensitivity to the 
factorization scale. 
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Work remains to make our joint resummation formalism a tool for phenomenological predic- 
tions in 1PI cross sections. From a numerical point of view, the resummed expressions involve 
extra integrals compared even to threshold resummation, so purely computational considera- 
tions make it more challenging to implement. Perhaps more significantly, it will be necessary to 
develop an appropriate matching formalism for large recoil. Nevertheless, we believe that the 
joint resummation formalism sheds valuable light on the reliability of perturbative calculations 
in hadronic scattering, and on the influence of nonperturbative effects. 
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A Status of the Refactorizations 



Specific arguments for the refactorized cross section, 
given in Ref. |34[ 



Eq. (|6]) for threshold resummation were 
>D was analyzed in some detail at one loop 



in Ref. || 



The Qr-refactored expression, Eq. 
although explicit arguments for its validity to all orders have not, to our knowledge, 
appeared in the literature. Eq. (0) is new to this paper. In fact, we believe that all three of 
these refactorized cross sections, Eqs. (||), (H) and @ are on a theoretical footing similar to that 
of collinear factorization, Eq. ([]]), but, of course, with different corrections as indicated above. 
We summarize below the ingredients of general factorization proofs for these relations, modeled 
after the arguments for collinear factorization given in Refs. Jl|, Q and summarized recently in 
Ref. p3j. The following reasoning applies explicitly to electroweak annihilation cross sections; 



the extension to single-particle inclusive cross sections is straightforward [50] 



Al. Leading regions and cut diagrams 

Each of the functions in the refactorized expressions for electroweak annihilation, Eqs. (||), (P) 
and (Q), corresponds to one of the subdiagrams in Fig. [TJ, which represents a general leading 
region, in the terminology of Ref. Q. A leading region in momentum space is one that gives rise 
to a contribution to the cross section that is leading power in the hard momentum Q. In general, 
leading regions can be classified by a set of on-shell virtual lines, whose vanishing denominators 
produce logarithmic corrections. In the following, an on-shell momentum is one whose invariant 
mass is much less than Q, and a soft momentum means one all of whose components are much 
less than Q in the specified leading region [Qj. 

Fig. [I] is in cut diagram notation, in which C represents a particular final state, with a 
graphical contribution to the amplitude on the left of C, and a contribution to the complex 



41 



conjugate amplitude on the right of C . In this notation, the cross section is a sum over all cuts 
C of forward-scattering diagrams G, consistent with the specified final state. In this case, all 
relevant final states contain an electroweak vector boson V, with momentum In Fig. [I], the 
subdiagrams Ja, Jb-, U and H include, respectively, on-shell lines with momenta collinear to pa 
(to be absorbed into V a / a , ip a /a or 1Z a / a ), lines with on-shell momenta parallel to pb (to Vb/b, i>b/b 
or TZb/b), lines with soft momenta (to U), and lines off-shell by order Q (to H). 

Following the notation of Ref . [jl| , we represent the contribution from region L of graph G to 
the cross section as 

dk\ dk~ B „ (c) /,+ ,-\ TT f d% f d% T (C) f u+ .N{w-Mn} 



c 

XI. 



^(tf.??),,.... >^(^) M - d") 



In this expression, we have suppressed the flavor labels (a and b in Eq. ([!]) and so on), and 
the corresponding Lorentz indices, that link the jets with the hard scattering. Eq. ( |114j ) is 
a representation of the most general leading-power contribution to the cross section da/d 4 Q, 
considered as an integral in the space of loop and final-state momenta. We shall argue that for 
each such leading region L, Gl may be rewritten as in Eq. @, @ or (pj), up to the corrections 
indicated in those expressions. For this discussion, it is convenient to introduce the lightlike 
vectors v M and u^, in the directions of the incoming momenta p a and Pb, according to 

u" = 4 = <W, ^ = — = V- (H5) 

Pa Pb 



Already at the first step in our factorization argument, Eq. ( |114j ), we have made the approxima- 



tion that the short-distance function if(/c^,/c^) depends on the large light-cone components of 
the quarks that annihilate to produce the electroweak vector boson. In particular, we identify 
[k'X/p^) = x a and (k^/Pb) — x b in Eq. (|1|). We know of no other approximation that can be 
extended beyond lowest order to formulate the short- distance function in a consistent fashion, as 
a collinear- finite quantity [TJ. Thus, we neglect transverse momenta in the calculation of H . This 
implies that in our factorized cross section, as in Eq. (|5|), the measured transverse momentum, 
Qt does not appear in the short- distance function directly, but only through its kinematic linkage 
( "recoil" ) with the total final-state transverse momenta of the jets J, and of the soft subdiagram 



U . Notice that the short distance function H in Eq. (|114|) need not be identical to the partonic 



hard scattering function a, which will in general absorb infrared-safe corrections from U. 
Al.l. The soft approximation and the eikonal function 

The desired factorized expression for each case discussed above reduces to an identity, if we make 

(c) (C) 

the following substitutions for the jet functions J\ and J B , 

4 c, fe, 9 i) < " n> - j^te.fe-^")'" 1 ^...^'...^. (lie) 



These substitutions are the "soft approximation" [pj. In the soft approximation, the momenta 
and polarizations of all soft gluons connected to J a are approximated by their minus components, 
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(a) 

(b) 



Figure 6: (a) Gluons attached to jet J a in the soft approximation, indicated by the dashed 
line. Color indices d and e are to be contracted with the hard scattering, (b) Result of the 
Ward identity discussed in the text, in which the jet is replaced by an eikonal (double) line with 
velocity vector, v M . Note the color trace on the remaining jet functions, normalized by d a , the 
dimension of the color representation of parton a. 

that is, the components moving opposite to the direction of lines in the jet. Similarly, for Jb, 
only plus components are kept. Notice that although the transverse components of soft momenta 
are neglected in the jet functions, no approximation is made in the soft function itself. 

Once the soft approximations ( |116| ) have been carried out, the decomposition of Gl, Eq. 



(|114j) , into the appropriate convolutions in transverse momentum and light-cone or energy frac- 



tion, requires only that we sum over all connections of gluons from the soft diagrams U to the 
corresponding jet subdiagrams. The graphical Ward identities of the theory ensure this result, 
illustrated by Fig. || ||. After the sum over diagrams, the coupling of the soft gluons to jets is 
replaced by their coupling to eikonal lines, as in Eq. flllD , in the directions of the jet momenta, 
which serve as sources for the gluons in the functions U. The Ward identities, which are essen- 
tially algebraic, also do not require an integral over soft-gluon momenta, which may be treated as 
fixed. The use of Ward identities is thus consistent with the restrictions on final-state momenta 
necessary to define cross sections at observed transverse momenta or near threshold. 

A1.2. Justification of the soft approximation 

We have now seen that a refactorization appropriate to each of our theorems is assured if the 
soft approximation can be justified. What must therefore be verified is the accuracy of the soft 
approximation, as it is used in each factorized or refactorized cross section. As discussed, for 
example, in Ref. the soft approximation fails, on a graph-by-graph basis, for J\' whenever 
Qe qe t T for any one of the soft lines that connect to U, and in J B C ^ whenever g+ <C q^ T for 

(C) 

one of the lines that connect J B with U. The central step in the proof of collinear factorization 
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Q is to show that all regions where the soft approximation fails cancel in the sum over final 

— (C) 

states. Equivalently, we must show that the q e integrals linked to J A , considered as contour 
integrals in the complex plane, are not pinched between coalescing singularities with qj <C qi^- 
For a fixed final state, however, it is easy to identify such pinches, between poles in the upper 
half-plane from the Feynman denominators of "spectator" lines, which carry large plus momenta 
into the final state, and poles in the lower half-plane from the denominators of "active" lines, 
which carry plus momenta into the electroweak annihilation [ 541, HI- The former may be thought 
of as characteristic of final-state interactions, and the latter of initial-state interactions. After 
the sum over final states, however, final-state interactions cancel, and all remaining singularities 
in qj from the function J a are in the lower half-plane. As a result, after the sum, the qj integrals 
are no longer pinched in the dangerous region, and the soft approximation may be carried out. 

For collinear factorization, Eq. ([!]), the cancellation of final state interactions is verified once 
the internal minus momentum integrals internal to jet J a, as well as the internal plus momentum 
integrals of the remainder of the diagram, are carried out. Details can be found in Ref. Q; 
here we shall only need the result that no other integrals are necessary. We must verify that 
corresponding arguments apply to the Qt, threshold and joint refactorized cross sections. 



A1.3. Qt refactorization 

An essential feature of the cancellation of final-state interactions in the collinear-factorized cross 
section is that it requires integrals over light-cone momenta only || . The cancellation is therefore 
consistent with fixed transverse momenta. The new feature of the Qt- refactorized Eq. ([5]), relative 
to the collinear-factorized cross section, Eq. ([[]), is simply that in the former the total transverse 
momentum of the soft and jet functions is frozen at —Qt- The arguments for factorization as in 
Eq. ([[]) are therefore adequate for refactorization as in Eq. ([5]), where Qt is the total transverse 
momentum carried by quanta that are absorbed by the short distance function \H\ 2 in Eq. ( |1 14| ) . 

Eq. (||) refers to a cross section in which the observed final-state momentum is that of a single 
vector boson only, with all real QCD radiation incorporated into the jets or soft function. In 
general, the short distance function, \H\ 2 , like the jet and soft functions, includes corrections 
associated with QCD radiation into the final state. The structure of the leading regions in 
Eq. ( |114[ ), however, ensures that these corrections are not singular when the total transverse 



momentum of the extra, "hard" QCD radiation vanishes. By construction, singular behavior is 
entirely in the jets and soft function. As a result, such regions, although of leading power in Q 2 , 



are not singular as a power at Qt = 0, when Qt is the vector boson momentum || pq| . These 
nonsingular contributions, which begin at next-to-leading order, are absorbed into the correction 
Y kt in the Qt- refactorized cross section, Eq. (0). 



A1.4. Threshold refactorization 

The situation for threshold resummation is a bit more subtle. In this case, we want to fix 
the energy of radiation from the jets and soft function in the hard-scattering center-of-mass 
frame. This clearly puts restrictions on the light-cone momentum integrals needed to ensure 
the cancellation of final state interactions. Near threshold, however, all radiation into the final 
state, including the radiation within the jets, is soft compared to Q, because near threshold 
the total energy of final state radiation is of order (1 — t)Q, with r = Q 2 /S. All radiation, 
including "spectator" gluons in the jets may thus be absorbed into the soft function U in Eq. 
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(|114j) . Corrections to this result are less singular by a power at r = 1 than the, leading, 1/(1 — r) 
behavior in perturbation theory. The jet subdiagrams that remain after this factorization are 
purely virtual, with only "active" lines, whose energies flow into the hard scattering [H]. In this 
case, as pointed out in Sec. 3.2 of Ref. ]34|], the same mechanism for the cancellation of final-state 
interactions in collinear factorization, Eq. (|1|), ensures that the soft function near threshold is 
free of overlapping soft-collinear logarithms, even when the light-cone integrals are restricted. 
Collinear singularities are present in the soft subdiagram, which is now the eikonal cross section 
discussed above, and they factorize in the usual way, as in Eq. fl37|) above (integrated over fcr)- 

We can now reorganize the cross section in moment space to derive Eq. (^]). We remove 
factorizing, purely collinear singularities from U by multiplying and dividing the eikonal cross 
section by eikonal jet functions ipf lk \ corresponding to the incoming partons, and defined by 
analogy to Eq. (|22|), in the appropriate (Laplace) transform space. The ratio of the soft function 
to eikonal jets in transform space, as in Eq. fl37D , is free of collinear singularities altogether. 
The eikonal jets combine with the virtual remainders from the original factorization to form the 
functions ip at fixed energy, near threshold, as in Eq. (|6]). 

We expect there to be singular but integrable remainders in the soft function, after the sum 
over final state interactions. These remainders are precisely the logarithms associated with the 
soft function, U |34j . Corrections resulting from the factorization of final-state radiation into U, 
which are O ([1 — r]°), are absorbed into Y th in Eq. (H), and begin at one loop. 



A 1.5. Joint refactorization 

The step from the threshold- ref actorized expression @, to the jointly refactorized Eq. @ is 
essentially equivalent to the step from collinear factorization to the Q-r-refactorized form (||). The 
extra ingredient is again to fix transverse momenta for the soft radiation and the jets, a condition 
that we have already argued is consistent with the cancellation of final-state interactions, and 
hence with the validity of the soft approximation. 

Of course, as emphasized in Refs. JI], [| , complete factorization proofs should rest on arguments 
based upon systematic subtraction procedures. This level of sophistication is still to be attained 
in hadronic scattering cross sections. Nevertheless, we consider the arguments outlined above to 
be adequate to justify an analysis based on the refactorization theorems in Sees. 2 and 4. 



B One-loop parton distributions and refactorization 

In this appendix we present for illustration, at one-loop accuracy, the various generalized dis- 
tributions of section 2.3. To keep the presentation brief, we will only present results for quark 
distributions. We begin with the familiar distributions, defined at fixed 'plus'-component of their 
momentum, and then turn to those at fixed energy. The latter, with transverse momentum fixed 
also, are then used to verify the joint refactorization (|9]) to one loop for electroweak annihilation. 
Illustrations of joint refactorization for the one-loop prompt photon cross section will be given 
elsewhere. 



45 



Bl. Quark distribution at fixed light-cone momentum fraction 



We give first the one-loop partonic light-cone distribution ( |10"D in n ■ A = A + = 0, n 2 = gauge. 
The familiar result is : 



<P q/q (x, /I, e) = 5(1 - x) - (i + ln(47r) - 7 *) ^ C t 



1 + ,3^ 



(117) 



B2. Quark distribution at fixed light-cone fraction and transverse 
momentum k 

Next we compute V q / q (x, k, e), Eq. (fL7l ), in n ■ A = gauge with n 2 ^ 0. For simplicity, we 
choose riy = 0. One finds 



V q/q (x,k,n,e) = 6(l-x)5 D - 2 (k) l + ^/' v (x,2p-n,e) 



1 — x 
k 2 



1 -e + 



2x z/ 



k 2 + (l-x)V 



2z/(l - a;) 



(k 2 + (1 - x)V) 2 j ' 



118) 



with ip q ) q V the same function of the gauge vector given below for the one-loop virtual correction 
to i/j q / q , and where v = (2p ■ n) 2 /\n 2 \. It is not difficult to check that 

d d - 2 kV q/q {x,k,fi,e) = _ j Q dk 2 (k 2 )^ V q/q {x,k,fi,e) 

= <f) q/q (x,ii,e) + finite, (119) 
i.e., that the singularities in V q / q match those of 0, up to finite (n- and /i-dependent) remainders. 

B3. Quark distribution at fixed energy 

For the one-loop energy distribution (0), the computation is described in detail in Ref. ||. We 
work in n ■ A = 0, = 5^° gauge and treat the l/(n ■ k) terms in the gluon propagator in 
principal value prescription. The result is 



ip q/q {x, 2p Q , e) = 5(1 - x) + ip { q / q R {x, 2p Q , e) + %}^{x, 2p Q , e) 
where p$ = Q/2, and the real and virtual contributions are given by 



(120) 



x Re 



a s (fi)C F Utt^Y T(2-e ) ( -1 
2tt V v j T(2-2 



eve 



1 + x — e(l — x) 
(l-x) 1+2c 



(121) 



7T 



5(1 -x) 



2?) (rr^l 111 -'iin + i'i 



- (sign(n 2 )) £ + | (sign(n 2 ) e 2 ^'' 



+ 4 (- + 1b(4tt) - 7£ 



(122) 



where again v = (2p ■ n) 2 /\n 2 \. ipQ is identical. The counterterm contribution to Vvq results 
from MS fermion wavefunction renormalization in the n ■ A = gauge, and the double pole in e 
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reflects an overlapping soft and collinear divergence. The expansion in e of (|120|) reads 

'a q Cp\ I 1 /l 



*P q/q (x,2p ,e) = 5(l-a0+(^) ]--(-+ ln(47r)- 7E 
+2 



1 + ^ 2 , 3 



ln(l — x) 




■ 1 " 




■ 1 ■ 


1 — X 


+ 


.1 — x. 


.1 — X. 



7T 



+ —5(l-x) , (123) 



plus non-singular terms. Note at O(e ) the appearance of the double- logarithmic plus-distribution 
[ln(l — x)/(l — + , which is a remnant of the 1/e 2 cancellation between real and virtual con- 
tributions. 



B4. Quark distribution at fixed energy and transverse momentum 

The one-loop result for the distribution lZ q / q (x, k, 2p , e) in Eq. ([0J) is 



TZ q/q (x, k, 2p , e) = 5(1 - x) 5 D - 2 (k) + K^x, k, 2p , e) + nfF(x, k, 2p , e) 



where 



n { ^ R {x,k,2p Q ,e) 

n^; v {x,k^Po,e) 

In these expressions, we define 
F(x, k, 2p , e) 



F{x, k,2p ,e)P^(x) 



k 2 (2p ) 2 (1 - x) 2 



«£CV /. 2 2V f -, _ 4k 2 
^ V 47r/i J ^ (l-x) 2 (2p ) 2 

\2 



-1/2 



2tt 2 

1 + x 2 — e(l - r) 



1 — x 



Tlf/q is identical to T\j/ q - Consistency requires that 



d d ~ 2 kn q f(x,k,2p ,e) 



TV 



1-e 



? '" '" X)2 dk 2 (k 2 )-^(x,k,2p ,e) 



r(l-e) 



(124) 

(125) 
(126) 



(127) 



128) 



which is straightforward to verify. Note that in TZ the double-logarithmic singularities in 1 — x 
and k are generated by the term 



P 2(l-X) 2 



-1/2 



2/i_2\-l-e 



dk\k 



1 — X 



1 - 



(l-x) 2 p 2 0/ 



(129) 



B5. Drell-Yan cross section 

The lowest order cross section for the reaction 

q(p 1 ) + q(p 2 ) -> 7*(?) 



(130) 
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reads in D = 4 — 2e dimensions 



da, 



qq 



Aircx 2 



dQ 2 d D - 2 Q T 3N c sQ< 



'.l-e)5 1 



Q 2 



6 D - 2 (QT) = aoS(l-z)S u --'(Q T 



D-2, 



The one-gluon radiative correction to this cross section, resulting from the reaction 

?(Pi) + ?(P2) -> 7* (<?) + ^(fc) , 



(131) 



(132) 



is straightforwardly computed, using for example the expressions given in [pqj . In terms of the 
quantities defined in (|127|) , the result is written compactly as 



<Pcrj£' R (z,Q 2 ,Q T ) 
dQ 2 d D ~ 2 Q T 



2a J r (z,Q T ,Q/y/z, e) 



22 



;i - 



(133) 



with P e given in ( |127| ). One can easily verify that upon integration over d D 2 Qt this expression 



gives the real-emission correction to the inclusive Drell-Yan cross section given in Eq. (88) of |BB 



B6. Joint refactorization at one loop 

Let us now illustrate the refactorization of the Drell-Yan cross section in Eq. @ at one loop. We 
will see that its singularities in Q T and 1 — z are accounted for by the distributions i = q, q, 
and the soft function U qq -. 

At one loop, the right hand side of Eq. (pj) expands into a sum of the one-loop expressions for 
its factors. Virtual corrections in Eq. (|9|) are exactly equivalent to those of Eq. (Q) in threshold 
resummation |J. Thus, we restrict ourselves to the real contributions, and check that the hard- 
scattering function found by expanding (§) to one loop: 

a ° k ™ - dQ 2 d°- 2 Q T {z ' Q Mt) 

-a \fcf/q R (z, Qt, Ql \[z, e) + T^f^ R (z, Qt, Q/ \fz, e) + U$' R (1 — z, Qt)] , (134) 

is free of singularities at Qt = and (after Qt integration) at 1 — z — 0. The soft function 
Uq^' R may be found from Eq. fl3?D. At one loop, it is quite simple to determine, because only 
the interference graphs contribute in the ratio, and even in these diagrams only the k^k u / (n ■ k) 2 
term in the gluon polarization tensor survives. The result is: 

Uf> R (1 - z, Q T ) = T{z, Q T , Q/^, e) Q2 ^_ . (135) 
Using Eqs. (|T2§, (J33J) and (PSD , we then find 

a = -*o Hz, Qt, Q/V~z, e) + % . (136) 



This is the desired behavior, because, when integrated over Qt, the result is nonsingular as z — > 1 
(see the discussion after Eq. (pj)). 
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